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Abstract 

An expansion in the number of spatial covariant derivatives is carried out to 
compute the ("-function regularized effective action of 2+1-dimensional fermions 
at finite temperature in an arbitrary non-Abelian background. The real and 
imaginary parts of the Euclidean effective action are computed up to terms 
which are ultraviolet finite. The expansion used preserves gauge and parity 
symmetries and the correct multivaluation under large gauge transformations 
as well as the correct parity anomaly are reproduced. The result is shown 
to correctly reproduce known limiting cases, such as massless fermions, zero 
temperature, and weak fields as well as exact results for some Abelian config- 
urations. Its connection with chiral symmetry is discussed. 
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I. INTRODUCTION 



Three-dimensional field theories find a direct physical application in the context of con- 
densed matter [|]] and also appear as the high temperature limit of 3+1-dimensional theories 
|| . From the theoretical point of view, odd- dimensional field theories have some remarkable 
properties. They are free from scale anomalies || and so the corresponding coupling con- 
stants do not run at one- loop @]. Gauge theories in odd-dimensional spaces naturally admit 
a local term which is of topological nature, namely, the Chern-Simons secondary character- 
istic class (CS). Such terms are invariant under topologically small gauge transformations, 
induce a parity breaking in the effective action and give a mass and spin to the gauge boson 
||. Moreover, the so called topological mass of the photon in QED3 is not renormalized 
beyond one loop order [MJ. In the non-Abelian case, the CS term may change by integer 
multiples of 2iri under large gauge transformations, implying a quantization condition on 
the mass to coupling constant ratio |J. The quantization of the CS coefficient is in com- 
plete analogy with the quantization of the topological Wess-Zumino-Witten action (WZW) 
in even- dimensional space-times [J|,f§. 

When fermions are coupled to the gauge fields, the CS term is induced by quantum 
fluctuations @-ll|- Preservation of gauge invariance requires then the presence of an appro- 
priate CS counter term which introduces a parity anomaly in the effective action ||. The 
parity anomaly has the same technical meaning as for axial anomalies in even dimensional 
theories, that is, the anomalous term breaks a symmetry that was present at the classical 
level. In the chiral case the fermions induce a WZW term and enforcement of vector gauge 
invariance introduces an axial anomaly Hl2| . The dimensional ladder relating chiral, parity 
and gauge anomalies for massless fermions in space-times of dimensions D+2, D+l and D 



has been elucidated in 



The relation between axial anomalies in D+l dimensions and 
effective actions in D dimensions using trace identities, both for zero and finite temperature 
and density, has been studied in [pT0| , |T6H T8|1 . It is also noteworthy that induced CS terms 



appear in four dimensional theories at finite temperature and density |19 

The exact parity-violating piece of the effective action for massless fermions in odd- 
dimensional space-times was obtained in []nj in terms of the 77-invariant. This result follows 
from the use of the (^-function regularization [|20] which automatically preserves strict gauge 
invariance under large and small gauge transformations. The //-invariant measures the spec- 
tral asymmetry of the corresponding Dirac operator, and through an application of the 
Atiyah-Patodi-Singer index theorem [21,14], it can be related to the CS term corrected by 
the index of the Dirac operator extended to one more dimension. Actually, the massless ex- 
act result holds for zero as well as for finite temperature. Indeed, it applies to fairly general 
space-time manifolds and, within the imaginary time formulation, the finite temperature 
amounts to a S 1 compactification in the time direction. 

For massive fermions there is no exact closed solution at zero temperature but some exact 
results are known, including the presence of a correctly quantized CS term (see e.g. p2],f2~3|j). 
At finite temperature, however, inconsistencies did appear in the first calculations of the in- 
duced CS term in 2+1 dimensions both in the Abelian and non-Abelian cases |T6||24| . Since 
the CS term is a local polynomial, the calculations to extract this term were usually done 
using a combination of perturbative and derivative expansions (see also [18| for a recent 
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calculation of this type at finite temperature and density, and |25| for perturbative calcula- 
tions without the low momentum approximation). The result of these calculations yields a 
CS term with a coefficient which is a smooth function of the temperature and this breaks 
invariance under large gauge transformations. The necessity of gauge invariance was empha- 
sized in ||26|| . In |?7| it was argued that perhaps the exact (i.e., non perturbative) coefficient 
is actually a stepwise function of the temperature so that only properly quantized values 
appear and in any case only those quantized values would contribute after the functional 
integration on the gauge fields (this observation, however, does not solve the problem when 
the gauge fields are external). Replacing the CS by its gauge invariant version, namely, the 
77- invariant, does not solve the problem either. This is because this quantity is gauge invari- 
ant at the price of displaying jumps of ±2 coming from its index contribution. However, if 
the spectrum of the Dirac operator is discrete and sufficiently well behaved, the renormalized 
fermionic determinant will be free from discontinuities under continuous deformations of the 
gauge configuration, therefore the exact effective action can only have jump discontinuities 
which are multiples of 2iri (or multiples of in if some eigenvalue vanishes during the defor- 
mation; this is the generic case for massless fermions). This leads again to a quantization of 
the coefficient in the 77-invariant. 

This apparent paradox has recently been solved: in Ref. PSI it is shown, in the exactly 



solvable 0+1-dimensional Abelian case, that the exact result is free of pathologies and con- 
sistent with gauge invariance. Large gauge transformations (which exists in this Abelian 
model due to the compactification introduced by the temperature) mix different orders of 
perturbation theory. In Refs. [^| the 2+1 dimensional case is considered. There it is shown 
that, using the (^-function prescription, the effective action at finite temperature can be 
regularized in a completely gauge invariant manner. Furthermore, the imaginary part of 
Euclidean the effective action is computed explicitly for some special configurations where 
separation of the time and space variables is possible and the gauge fields are effectively 
Abelian. In Ref. |5D| another calculation is carried out for such configurations, this time by 
integration of the chiral anomaly of the associated two dimensional problem (the calculation 



is extended to arbitrary odd dimensions in |31j|). These results check invariance under both 



large and small gauge transformations and display the expected parity anomaly. Related 
work, both on odd and even dimensional space-times, prompted by the insights in those 
papers can be found in [|32|-|3"4|j . 

In the present work we study the Euclidean ^-function regularized effective action of 
fermions in 2+1 dimensions in the presence of arbitrary external non- Abelian gauge field 
configurations and at finite temperature. The calculation is done not only for the imaginary 
part, which contains the topological terms, but also for the real part. 

Besides the ^-function regularization, the main ingredients of the calculation are, first, 
the use of a Wigner transformation to represent the operators and second, an expansion 
in the number of spatial covariant derivatives. The Wigner representation method was 
introduced in [^| and is closely related to the method of symbols for pseudo-differential 
operators. (An improved symbols method and many references can be found in ||35|| .) The 
method is adapted here to treat the finite temperature case, and has been used also for 1+1- 
and 3+1-dimensional fermions at finite temperature in p3 and for non-local Dirac operators 
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in 1 36|| . The combination of (-function and Wigner transformation yields a computational 
setup which is suitable to obtain the effective action using different expansions. [] 

The other ingredient is the use of a derivative expansion. In order to preserve gauge 
invariance, it is, of course, necessary to expand in terms of the covariant derivatives. This 
would be sufficient at zero temperature, however, at finite temperature because the frequen- 
cies take on discrete values only, an expansion in the number of time-like covariant deriva- 
tives breaks invariance under large gauge transformations and so such derivatives have to be 
treated non-perturbatively. Therefore, we will consider an expansion in the number of spatial 
covariant derivatives only. The advantages of such an expansion are, first, that it preserves 
gauge and parity invariances, i.e., terms of different order do not mix under such transfor- 
mations, and second, each order increases the degree of ultraviolet convergence. Therefore, 
computing the terms which are ultraviolet divergent, as we will do, yields all contributions 
which may contain anomalies and multivaluation under large gauge transformations. In 
particular, topological terms are included there. The remaining terms are strictly gauge and 
parity invariant. 

For technical reasons, the calculation is carried out in the gauge doAo = which is shown 
to exist for any given configuration, although it is not unique. We check that different 
d A = gauges yield compatible result. The calculation is not based on integrating the 
axial anomaly in two dimensions or similar ideas. Except for the use a particular gauge, it is 
a direct and systematic calculation which uses only the definitions of (-function, derivative 
expansion, etc. This allows to obtain both the topological and non-topological terms. The 
expansion can be carried out to any order in principle, although, of course, higher orders 
require an increasing amount of work which quickly becomes prohibitive. As mentioned, in 
the present work we compute the terms with ultraviolet divergences. This amounts to zero 
and two derivatives for the real part and two derivatives for the imaginary part. We check 
that our result is consistent with gauge invariance, more precisely, it transforms as the CS 
term, and yields the correct parity anomaly, which is shown to be temperature independent. 
Further, it reduces to the result in |30| for the same configurations, has the correct zero 
temperature limit and reproduces the exact result for massless fermions. 

Another related expansion is also considered which consists in carrying out a further 
expansion in powers of the time-like covariant derivatives which are inside commutators (as 
will clear below, the space-like covariant derivatives were already inside commutators due 
to gauge invariance since the space is not compactified) . This new expansion still preserves 
gauge and parity invariances and higher order terms are increasingly convergent. At leading 
order, the result is a simpler expression which still contains all topological terms, saturates 
the parity anomaly and possesses the same multivaluation of the full effective action. This 



It should be noted, however, that the ("-function is not particularly convenient to compute the 
effective action in odd-dimensional space-times at zero temperature. This is because there is no 
pathology free way to choose the branch cut in z s so as to avoid both the spectrum of the Dirac 
operator, which is continuous under most expansions, and the spectrum of the parity transformed 
Dirac operator. This regularization is suitable at finite temperature since then the spectrum 
becomes discrete for practical purposes. 
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result also reproduces that of 3(J , as well as the zero temperature and zero mass limits. The 
same expansion has been used in [34| for even-dimensional fermions at finite temperature. 

The paper is organized as follows. In section [TT] some general considerations are made re- 
garding the effective action, its relevant symmetries and the Wigner transformation method. 
In section [TT1| the 0+1- dimensional case is revised. This allows to illustrate some exact re- 



sults as well as the method before going to the three-dimensional case. In section |V| the 
main results of the three-dimensional problem are presented. For the leading term of the 
imaginary part, which is the finite temperature version of the induced CS term, we check 
that it possesses all expected correct properties under symmetry transformations and limit- 
ing cases. In section [V] we present the explicit computation of the effective action at lowest 
order, that is, without derivatives. This term is purely real. Section [Vl| contains the compu- 
tation of the terms which are real and with two derivatives. Section |V1J] contains the similar 
calculation for the imaginary part. Appendix [A] proves the existence of the A -stationary 
gauge and finds the allowed gauge transformations within that gauge. Finally, Appendix [B| 
introduces some mathematical results used in the main text. 



II. GENERAL CONSIDERATIONS 

A. The Dirac operator 

The most general Dirac operator in a flat three-dimensional Euclidean space has the 
form 

D = llL D li + M (1) 

where = + A M is the covariant derivative and A^x) and M(x) are matrices in some 
internal space to be referred to as flavor. Our convention for the Dirac matrices will be 

7j = 7^, {7 M , 7j = 2<V- ( 2 ) 

In 2+1 dimensions there are two inequivalent irreducible representations of the Dirac algebra 
labeled by r\ = ±1, defined as 7o7i72 = if}- For instance, 7 = r]T 3 , 7 12 = r 12 where T; 
stand for the Pauli matrices. To ensure unitarity in Minkowski space, is required to be 
ant i- Hermit ian and M Hermitian. At finite temperature T, the Dirac operator D acts on 
the space of fermionic wave functions ip(x), which are antiperiodic in the Euclidean time 
direction with period j3 = 1/T. Correspondingly, the external bosonic fields A^ and M are 
periodic. Further restrictions on D comes by imposing continuity in all involved functions 
(fermionic wave functions, field configurations, gauge transformations, etc). In addition, 
we will occasionally assume that the space is compactified so that it is topologically a two 
dimensional sphere. This implies a topology M3 = S 2 x S 1 for the space-time. The class of 
Dirac operators to be considered in each case must be a subset of the above and should be 
sufficiently large as to be invariant under the relevant symmetries of the problem. 



B. The effective action 

The unrenormalized partition functional is 
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Z(D)= I VipVip e~ J d3x ^ 



(3) 



Formal integration of the fermions gives Z = Det(D) and thus the unrenormalized effective 
action is 

^w(D) = - logZ(D) = -Tr log(D) . (4) 

Upon regularization and renormalization of the ultraviolet divergences, a well-defined and 
finite effective action W is obtained. However, the renormalization procedure is not unique 
and thus there will be several functions W(D) all of them qualifying as valid effective ac- 
tions for the same original action / d 3 xi/jDip. They are constrained to reproduce the same 
ultraviolet finite terms since such terms can be computed without regularization and are 
therefore unambiguous. This implies that the fourth order variation should be common to 
all VT(D). This is because for arbitrary commuting variations of D 

8 1 5 2 5 3 8 4 W(B) = Tr^DD^DD-^DD-^DD- 1 ) + 5 permutations (5) 

is ultraviolet finite in 2+1 dimensions. As a consequence, the (renormalized) ultraviolet 
divergent terms, responsible of the ambiguity in the renormalization, must vanish after a 
fourth order variation. Thus the allowed ambiguity introduced by the renormalization is 
a local polynomial action of canonical dimension at most three constructed with the fields 
M, and their derivatives. This is also clear in perturbation theory since all Feynman 
diagrams become convergent after four derivatives in the fields or the external momenta.^ Of 
course, in order to compare the effective action for different Dirac operators, it is essential to 
use the same renormalization for the whole class of Dirac operators D considered^] Only in 
this case can be meaningful the statement that W(D), as a function of D, is unique modulo 
a local polynomial. All these remarks are well-known but they will be needed later and non 
trivial consequences will be extracted from them. 

A set of renormalizations are based on the spectrum of the Dirac operator. Let Dip n = 
\ n ip n be the eigenvalue equation for D. The ultraviolet divergent part of the spectrum comes 
from the term 7^*9^ of the Dirac operator and thus it lies on ±ioo. In practice we will take 
M(x) — m, a constant c-number, and in this case the spectrum will lie on the line m + iR. 
The bare effective action can be expressed as 

W / bare = -E 1 °g(An). (6) 

n 



2 There is a subtlety at finite temperature. In this case the frequency takes discrete values only and 
this prevents to take a derivative with respect to it. Nevertheless, since the ultraviolet convergence 
also increases by taking finite differences, the argument applies also for finite temperature and the 
allowed ambiguity is still a local polynomial. 

3 As an extreme case of ignoring this restriction, consider adding to the effective action a different 
"constant" for each D. This would result in a completely arbitrary function W(D). 
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In this form, the sum over eigenvalues is ultraviolet divergent and must be renormalized. 
However, the renormalization only affects the ultraviolet part of the spectrum so any par- 
ticular eigenvalue still contributes with — log A n to the renormalized effective action. This 
shows that qualitatively W(D) is to be understood as a many-valuated function defined 
modulo 27rz on the manifold of (sufficiently regular) Dirac operators to be denoted by Ai. 
The branching points correspond to the Dirac operators which are singular, that is, with a 
zero eigenvalue. 

A mathematically well founded choice for the renormalization prescription, which we will 
adopt here, is the (-function version of W, namely, 

iy(D) = -£A>gA 

n 

(The same branch is to be taken for the two functions z s and log(z) and all the eigenvalues.) 
This expression is to be understood as an analytical extension on the variable s from Re(s) < 
—3 and a sufficient condition to be well defined is that the eigenvalues are all on the same 
Riemann sheet of the logarithm function and off the branch cut |[20|| . Once the branch cut, T, 
for the function — z s \og(z) on the complex plane is chosen, VT(D) becomes completely well 
defined but it will, of course, have a spurious discontinuity along the cut. More precisely, 
the cut on the complex plane will produce on Ai a branch cut stemming from each of the 
singular Dirac operators. The corresponding cut manifold will be denoted by Air- On the 
complex plane and upon analytical extension, two different choices of the cut, I\ and T2, will 
yield the same function —z s \og(z). This is no necessarily the case on Ai. The two choices 
of the cut will yield the same function W(D) upon analytical extension if and only if r\ and 
T 2 can be deformed into each other crossing at most a finite number of eigenvalues. Because 
the ultraviolet sector of the spectrum lies on ±ioo there are essentially only two inequivalent 
prescriptions, namely, choosing the branch cut of the logarithm along the negative real axis 
or along the positive real axis. The two choices will be labeled by a = +1 and a = — 1 and 
the two cut manifolds by Ai a= \ and Ai a= -i, respectively. 

Analogously to the theory of many-valued functions on the complex plane, an alternative 
way to achieve both one-valuedness and analyticity of W(D) is to work on the simply 
connected Riemann surface manifold associated to Ai (i.e., its the universal covering space) 
which will be denoted by Ai. From the previous discussion it follows that the two functions 
W a =±i(D) defined on their common Riemann surface Ai by analytical extension represent 
two different ((-function) renormalizations of the effective action and thus must differ by a 
local polynomial function on Ai. 

Another comment refers to symmetries. As usual, a symmetry is defined as any trans- 
formation of the Dirac operator which leaves invariant the class of allowed Dirac operators 
and which can be compensated by a corresponding transformation in the fermionic wave 
functions so that the action / d 3 xipDif) remains unchanged. According to this definition, 
symmetries leave invariant the classical effective action (defined as the action evaluated at 
its minimum with respect to the fermionic fields). A non- vanishing variation of the quan- 
tum effective action implies a quantum anomaly for the symmetry. The anomaly must be 
a local polynomial. Indeed, if W(D) is any definition of the effective action and D n is the 
transformed Dirac operator with Q independent of D, the function W(D Q ) also qualifies as 
an admissible definition of the effective action of D. Therefore, the anomaly iy(D n ) — W(D) 
must be a local polynomial. Instances of this are the chiral anomaly in even dimensional 



s=0 



(7) 



7 



space-times or the parity anomaly in odd dimensions. It also applies to the 2irin multival- 
uation under large gauge transformations.^ 

In practice, the ^-function renormalization is carried out by introducing the quantity 

n s (D)=Tr(D s )=Y,K, (8) 

n 

which is ultraviolet finite if Re(s) < —3 and a meromorphic function of s with simple poles 
in s = —3, —2, —1 pD| . Then, the (^-function prescription is 

W(D) = — n.(D) _ n . (9) 



ds 



s=0 



The precise definition of the function z s requires cutting the complex plane along a ray 
characterized by an angle 9. Applying Cauchy's theorem [^0| 

fi sr (D) = -Tr / . (10) 

The trace Tr is to be taken in the fermionic Hilbert space which is the tensor product 
of space, time, flavor and Dirac. The integration path T follows the ray 9 starting from 
infinity, encircles zero clockwise and goes back to infinity along the ray 9 — 2tc, so that if the 
spectrum of the operator D were bounded T could be deformed to enclose it anticlockwise. 
The choices 9 = n and 2n correspond to a = +1 and —1 respectively. As noted, any value 
7r/2 < 9 < 37r/2 gives the same function Q S)Cr=+ i on the Riemann, surface, and similarly 
any 37r/2 < 9 < 5n/2 yields the same function f2 Sj(J= _i. Besides, The same effective action 
follows from taking 9 or 9 + 2™ for integer n. (See Appendix |B|) 



C. Gauge transformations 

Gauge transformations are defined as D u = f/ -1 D?7, where U(x) is a matrix valued 
function acting as a multiplicative operator on the fermionic wave functions. More explicitly, 
D u = ry li (d /t + A%) + M u , with 

A"(x) = U- 1 (x)(d fl + A^x^Uix) , M u (x) = U- 1 (x)M(x)U(x) . (11) 

U(x) belongs to some gauge group G which is a subgroup of XJ(Nf), Nf being the number 
of flavors. Correspondingly must be in the Lie algebra of G and likewise the class of 
matrices M(x) must also be closed under gauge transformations. (This is trivially satisfied 
if M is just a c- number.) Besides, U(x) must be continuous on the space-time manifold and 
in particular periodic function of x . 



It follows that it is always possible to define effective actions free from anomalies associated to 
compact groups, namely, by taking an average over all symmetry transformed configurations. This 
applies in particular to parity transformations, however, the average of the two effective actions 
defined modulo 2ni will have an ambiguity itt in general, so gauge invariance is not ensured. 



8 



Being similarity transformations, gauge transformations leave the spectrum of D invari- 
ant. As a consequence, the functions W a= ±i{D) defined on the cut manifolds M. a -±\ are 
strictly gauge invariant by construction since these functions depend solely on the spectrum 
of D. In this sense, the (^-function renormalization prescription trivially preserves gauge 
invariance. On the Riemann surface M. the statement is less trivial. Let D be continuously 
deformed along a path connecting two field configurations which are related by a gauge 
transformation. This will induce a corresponding bounded path on the complex plane for 
each eigenvalue and, since the spectrum is unchanged at the end, the net effect on the eigen- 
values can be at most a permutation. During its walk, a finite number of eigenvalues will 
cross the branch cut, each time adding a net ±27r«, to the effective action.^ In summary, 
the functions W a =±\{D) extended to M. are gauge invariant modulo 2ni. It should be noted 
that, by continuity, the effect of a gauge transformation may only depend on the homotopy 
class of the gauge transformation and on the homotopy class of the initial field configuration. 
In particular, only topologically large gauge transformations may change the effective action 
by a non vanishing integer multiple of 2iri. We are assuming that no branching point lies 
on the path followed by the eigenvalues. In the massless case this cannot be avoided since 
the spectrum is purely imaginary. In this particular case the variation in W a= ±\{D) is a 
multiple of in instead. 

For convenience, the calculation of the effective action will be carried out in an A - 
stationary gauge, i.e., a gauge such that d A (x) = 0. Such a gauge always exits for any 
given gauge configuration, therefore there are no loss of generality by making this choice. 
This result is proven in Appendix [A|. In the same appendix, it is shown that the gauge 
transformation needed to bring a gauge field configuration to an A -stationary gauge can 
always be chosen to be topologically small in the Abelian case. On the other hand, for 
simply connected gauge groups large transformations may be necessary, depending on the 
given gauge configuration. 

To check gauge invariance of our calculation, it will also be necessary to find the most 
general gauge transformations which preserve the gauge fixing condition. As shown in 
Appendix |A|, these transformations have the form 



Here Uq(x) is an arbitrary time-independent gauge transformation. A is a time-independent 
function taking values in the Lie algebra of the gauge group. It is only restricted by continuity 
of U(x), which requires 



Therefore, the most general gauge transformation preserving the Ao-stationarity condition is 
composed of a particular time- dependent transformation subject to the condition in eq. ([13]) 
followed by a stationary gauge transformation. The corresponding zeroth component of the 
transformed field takes the form 



Recall that each particular eigenvalue adds — log A n to the renormalized effective action; the 
renormalization, lim s _»oX) n A^logA n , leaves no trace on any given bounded region of the complex 
plane. 



U(x) = exp(-x A (x)) exp(x (A (x) + A(x)))U (x) . 



(12) 




(13) 
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A%{x) = U 1 (x)(A (x) + A(x))U (x) . (14) 

The time- dependent transformations e - xoA °e x °( Ao+A ) will be called discrete gauge transfor- 
mations for the stationary field Aq (since generically they form a discrete set due to eq. fljjD ). 

In general the condition on A can be simplified. From eq. ([13]), it follows that both Aq(x) 
and A(x) must commute with exp(/3A (x)). If the spectrum of exp(/3A (x)) is nowhere 
degenerated, this unitary matrix takes a diagonal form in an a?-dependent basis in flavor 
space that is essentially unique, thus A (x) and A(jc) must also be diagonal in the same 
basis and therefore they commute with each other. In this case the condition on A becomes 

[A (x),A(x)]=0, exp(/3A(aj)) = 1 . (15) 

The second condition means that the eigenvalues of A(aj) are of the form Xj = 2irinj//3, for 
integer rij. Such integers are cc-independent by continuity. Correspondingly, in this case, 
the allowed discrete transformations take the simpler form 

U(x) = exp(x A(x)) . (16) 

This form is the generic one. It applies if the spectrum of exp(j3Ao(x)) is nowhere degener- 
ated, or if the points of degeneracy can be resolved by continuity. 

The stationary gauge transformations are always topologically small since the second 
homotopy group is trivial for any Lie group. Whether the discrete transformations are 
topologically large or not depends on several factors. If the gauge group is not simply 
connected, for instance the Abelian case, U(l), and A is not zero, such transformations are 
necessarily large, since the loop around S* 1 for each x cannot be contracted to the identity. 
If the gauge group is simply connected, for instance SU (Nf), the discrete transformation 
may be large or small depending on A(x). For instance, if A is ^-independent, or becomes 
so after a continuous deformation, the discrete transformation will be topologically small. 
In particular, this holds if A(x) can be diagonalized using a similarity transformation which 
is continuous on S 2 , since a diagonal A must be ^-independent due to eq. (|15|). On the other 



hand, to see that there are large discrete transformation in general, consider the discrete 
transformations from S 2 xS 1 into SU(2), 

U(x,Xq) = exp ^™ X ° irx^J , neZ (17) 

where r are the Pauli matrices and x belongs to the unit sphere S 2 in R 3 . U(x) covers SU(2) 
2n times (it is 2n to one and it maps a positively oriented volume element on S 2 xS 1 into a 
positively oriented volume element on SU(2)). Equivalently, one can compute the normalized 
Wess-Zumino-Witten 3-form (relevant in chiral theories in 1+1 dimensions) which is closed 
and thus invariant under deformations, 

r wzw (C/) = I ti{R 3 ) = -Ainn , R = U^dU. (18) 

127T JM 3 

Therefore, such discrete transformations are all homotopically inequivalent and large for non 
vanishing n. The Wess-Zumino-Witten action appears naturally in this context since the 
Chern-Simons action 
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W CS (A) = -^[ tr(AF-\A 3 ), (19) 

47T JM 3 i 

(where F = dA + A 2 ) evaluated for a pure gauge field A = XJ~ X &U = R is just r wzw ([7). 
As noted in the introduction, the Chern-Simons action can only change under gauge trans- 
formations which are topologically large. 

D. Parity transformation 

The parity transformed of D will be defined as D p = 7^(9^ + A P ) + M p with 

A p {x) = -A {x p ) , A(x) = +A(x p ) , M p {x) = -M(x p ) , x p = {-x , x) . (20) 

This is the definition taken in ||14|| . All definitions of parity consist in reversing an odd number 
of coordinates and are related by a similarity transformation thereby being equivalent. 

The fermionic field transforms as ip p (x) = r jQip(x p ) and the eigenvalue equation becomes 
D p ip p = —\ n %l) p . Since the eigenvalues are changed, W(D) needs not be parity invariant 
within the ^-function regularization. This allows for a parity anomaly which is defined as 

A P (D) = W(D) - W(D P ) . (21) 

As noted at the end of subsection [II B| , the parity anomaly is due to the ultraviolet diver- 
gences present at the quantum level and hence is a local polynomial in the fields. The same 
result follows from noting that the two choices for the cut, o = ±1, are related by a parity 
transformation and thus 

W a (D p ) = W^(D) . (22) 

As a consequence, the anomaly can also be written as Ap(D) = W a (I)) — W_ CT (D), which, 
as noted above, is a local polynomial action on M.. 

Under parity, the action can be split into a parity preserving part plus a parity violating 
part, or equivalently, a a-even plus a cx-odd component 

W even (D) = i(W ff (D) + W_ ff (D)) , WWD) = \Ap(D) . (23) 

Since W dd is just a local polynomial, it can be removed by counterterms to end up with a 
parity preserving effective action. However, this can be incompatible with gauge invariance. 
By construction W a and Ap are gauge invariant modulo 2iri, but M^even and W Q dd can 
separately change by a multiple of in in general. When this happens, adding counterterms 
to remove W Q dd spoils gauge invariance of the effective action. Parity and gauge invariances 
cannot be both enforced in general 0. 

Because the parity anomaly is due to divergences in the ultraviolet sector of the theory, 
it is to be expected that it is essentially unchanged by introducing a finite temperature, 
which only affects the infrared sector through the periodic boundary conditions. The same 
argument applies for other quantum anomalies. In particular, it is well established in the 
literature (e.g. |37| ) that the axial anomaly takes the Bardeen form independently of the 
temperature. Remarkably, the analogous result can be proven for the parity anomaly, even 
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without a detailed calculation, namely, the form of the essential parity anomaly at finite 
temperature is given by the Chern-Simons action. By essential it is meant the part of the 
anomaly which is not removable by gauge invariant counter terms (i.e., the analogous of the 
Bardeen anomaly in the chiral case). To prove this result, we start with the parity anomaly 
at zero temperature. It is given by the Chern-Simons action which, of course, is a local 
polynomial, parity odd and gauge invariant modulo 2Tri. Therefore, at zero temperature, 
W dd = ±^W^cs- At finite temperature, W dd will be given by ±^Wcs plus a local polynomial 
with a smooth temperature dependence. By continuity this local polynomial must be strictly 
gauge invariant (not just modulo iir) and thus it can be removed by a counterterm without 
destroying gauge invariance, leaving Ap = iW^cs & t an Y temperature. This conclusion is 
confirmed by our detailed calculation below and also follows from eq. (5) of |H| (noting the 
different definitions of parity violating terms in both works) . The same argument applies to 
any other smooth dependence such as that introduced by a finite density or a change in the 
mass of the fermion. 



E. Pseudo-parity transformation 

It is also of interest to introduce a pseudo-parity transformation M, defined as 

Aq{x) -> -A (x p ), A(x)^+A(x p ), M(x) -> +M(x p ) . (24) 

It is clear that the pseudo-parity odd component of the effective action is that with an 
odd number of zeroth Lorentz indices and thus is the component containing the Levi-Civita 
pseudo-tensor. It also corresponds to the imaginary part of the Euclidean effective action. 
The CS action is pseudo-parity odd. The (^-function prescription yields a parity anomaly 
both in the pseudo-parity even and odd components, but only the latter is essential in 
the sense that it cannot be removed without breaking gauge invariance. This is because the 
even component can be renormalized preserving gauge and parity invariances simultaneously 
13|)14j|. Note that in the massless case, parity and pseudo-parity transformations coincide. 



F. Wigner transformation method 



The (asymmetric) Wigner transformation method is described in great detail in [23 
The corresponding formula for the zero-temperature (i.e., uncompactified) case is given in 
that reference and only a variation is needed to cope with the compactification in the zeroth 
direction. We concentrate on the Hilbert spaces of periodic and antiperiodic wave functions 
of time defined on S 1 . Let us introduce the fermionic and bosonic Matsubara frequencies 

2ni 1 lixin 
u n = —{n + -), Kn = ^~- ( 25 ) 

The states \u) = e ux ° are periodic if u = K n and antiperiodic if u — u n . In both cases they 
are normalized to y^fi. A general operator acting on the antiperiodic space will be of the 
form f(x , do) where x and d denote the time and frequency operators respectively, and / 
is a periodic function of x . Then the trace of / is given by 
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Tr(f(x ,d )) = ff^(w n \f(x ,do)\wn) 

P nGZ 

= ^E(°l e ^°/(^,4)e^ |0) 

P ngZ 

= 4E( l/(^o,4 + ^n)|0). (26) 

P ngZ 

Here |0) is the zero frequency state. Note that |0) is periodic rather than antiperiodic. For 
the non compact coordinates the treatment is entirely similar, Di is replaced by D{ + pi, 
there is an integral over momenta and |0) is the zero momentum state |23[| . Therefore, using 
eq. (|T0D, the ^-function can then be written as 

a -^-l§?^L^/^T^ 0) ' (27) 

Here p = ik (where k 6 R 2 ), po = u n and |0) represents the spinless and flavorless zero 
energy-momentum state normalized as (x|0) = 1. In particular <9^|0) = (0|<9^ = 0, and also 
(0 1 0) = / d 3 x, and moreover, if f(x) is a multiplicative operator, i.e., it does not contain 
derivatives, 

(0|/(ac)|0) = Jd 3 xf(x). (28) 

The formula just given for f2 S; r(D) is completely finite in the ultraviolet sector and exact 
if a space-time topology F^xS 1 is assumed. 

G. Covariant derivative expansion 

We will find convenient to use an A -stationary gauge to compute the effective action. 
This is allowed due to gauge invariance, however, since the A -stationary gauge is not 
unique, we will have to check that the result is the same for any of the gauge copies, i.e., 
invariant under time-independent and discrete gauge transformations. It should be noted 
that the integrand in eq. (^) is not directly gauge invariant due to the presence of |0); gauge 
invariance requires integration over momentum and sum over frequencies ||23|| . Roughly 
speaking, invariance under time independent transformations means that Di appears only 
inside commutators and this is equivalent to say that the result is unchanged if Di is replaced 
by Di + cij, a being an arbitrary ^-independent c- number, but precisely this invariance is 
ensured by the integration over momentum since a can always be absorbed in p. Similarly, 
invariance under the discrete gauge transformations, A — > A + A, is ensured by the sum 
over discrete frequencies, since A can be absorbed in uo n (recall that the spectrum of A 
is quantized). We note that a recently proposed technique, due to Pletnev and Banin 
p5| , allows to bring the expression to one where Di appears inside commutators from the 
beginning, i.e., prior to momentum integration. Presently this technique applies to covariant 
derivatives in non compact directions. 

In order to carry out the traces, sums and integrations implied in eq. (p7|) we will perform 
a formal expansion in powers of the spatial covariant derivatives, Di, i = 1,2. However, it 
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is important to remark that such an expansion is not only formal; it is meaningful for the 
effective action functional itself beyond the particular computational procedure followed. 
It actually means to change the bosonic field configuration to A (Xx, x ), XA(Xx,x ) and 
M(Xx,xq), and then count powers of A. Since this scaling preserves full gauge invari- 
ance (namely, by scaling U too) the gradient expansion is gauge invariant order by order. 
Formally, the scaling corresponds to the replacement Di — ► XDi and gauge invariance is 
preserved since Aa« can still be compensated by a change of variables in fcj. The situation 
is completely different for a formal expansion in powers of D Q . First of all, it would only be 
formal, since scaling xq —>■ Xxq violates the boundary conditions for most field configurations 
(because (3 is not scaled). At a formal level, scaling D would imply a similar scaling in A 
which will no longer be properly quantized and it could not be compensated by a shift in 
to n . 

Due to rotational invariance, the proposed gradient expansion contains even orders only. 
In this work we explicitly work out the zeroth and second order terms, Wo(D) and W2 (D) 
respectively. Since all higher order terms are ultraviolet finite and the different orders do not 
mix under gauge or parity transformations, all anomalies are contained in these two terms. 
The pseudo-parity odd component starts at second order since the Levi-Civita pseudo-tensor 
requires the presence of at least two gradients. Thus the leading pseudo-parity odd term, 
W^~(D), is responsible for all essential anomalies in the effective action, whereas W and 
may contain at most removable anomalies. 

There is another expansion which is also of great interest in the present context. As will 
be clear from the calculation, the final result at each order in the spatial gradient expansion 
comes as a function of A and V which is defined as [Do, ]. (Note that A and V commute 
in a A -stationary gauge.) It is the expansion in powers of A the one that would break 
gauge invariance under discrete gauge transformations, and so A should be treated non 
perturbatively. On the other hand, expanding in powers of T>q does not break any symmetry 
of the problem (excepting, of course, Lorentz invariance in the limit of zero temperature). 
Each extra power of T>q or Di increases the degree of convergence in the ultraviolet sector. 
Therefore, at finite temperature, it is natural to consider a double expansion in powers of D^ 
and T>q which play a similar role as a standard space-time gradient expansion which is often 
used at zero temperature. (In fact, Di always appears inside of commutators and so the 
expansion in powers of Di can be seen also as an expansion in powers of T>i = [Di, ]). This 
same idea has been applied in the case of 3+1-dimensional fermions at finite temperature 

III. THE 0+1-DIMENSIONAL CASE 

Before going to 2+1 dimensions, we will study the simpler 0+1-dimensional case. This 
will allow us to illustrate some of the previous remarks as well as the method. 

In 0+1 dimensions, the Dirac operator is 7o-D + m with 70 = 77 = ±1, a 1 x 1 dimensional 
matrix in Dirac space. We have already taken M(x) — m, a constant c-number. We start 
by fixing the gauge to be Ao-stationary. For convenience we will use units f3 — 2. A direct 
application of the Wigner transformation formula (i.e., eq. (E7|) adapted to 0+1 dimensions) 
gives: 
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1 r dz 1 
a r (D) = — Y / z s tr(0| — — 10) 

= _W— z s tr( - ^ (29) 

^7r2vri ^(w n + A ) + m - z) 

We have replaced D by v4 since everything is time independent and <9 vanishes on |0), 
and (0|0) = (3 = 2 has been used. 

In order to sum over frequencies, we first transform the series in f2 S) r(D) into a convergent 
one by performing a subtraction 

n sT m) = - I —z s tvY ( - A ; ) • (so) 

Jr2m ^ \r](u n + A ) + m - z r)(u n + A ) + m J 
This is justified due to the identity 

X^«W = °' (31) 

for any polynomial q(z), which comes from closing T with the circumference at infinity. 
This and similar formulas are to be understood upon analytical extension from sufficiently 
negative Re(s). Then, the following formula can be used 

n \Xi + U n X 2 + UJ n / 

and the sum over frequencies yields 

r 

VL S r (D) = - / — -z s tr tanh(m + r]A - z) . (33) 

The final step is to carry out the integration over z and apply —d/ds\ s= o to obtain the 
effective action. This can be done using the functions 

r dz 

tt r (tu, s) = - —z s tanh(uj - z) . (34) 
Jr 2m 

The properties of these and related functions are summarized in Appendix ||. As shown 
there 

n a {uj,o) = 

fi>,0) = log(e- 2 ^ + l). (35) 

Here the prime refers to derivative with respect to s and a = ±1 refers to the two possible 
inequivalent choices of T. Two remarks can be made regarding the first of these equations. 
First, the vanishing of f2 CT (u;, 0) guarantees that taking T along the ray 9 or 9 + 27m yields 
the same final result for the effective action. Second, it also implies that there is no scale 
anomaly in 0+1 dimensions. Indeed, the scale anomaly can be associated to a non trivial 
dependence of the effective action on some mass parameter M introduced by using (z/M Q ) s 
instead of z s in the definition of the ^-function. Such dependence cancels in odd dimensions 
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The final exact result for the effective action in 0+1 dimensions, restoring arbitrary units, 
is thus 



W a (D) = -tr log(e-^( m+r?j4 °) + 1) . (36) 

Recalling that under gauge transformations the eigenvalues of A may change at most by 
integer multiples of 2ni//3 it follows that Wo-(D) is manifestly gauge invariant modulo 2ni on 
the Riemann surface M. and strictly gauge invariant on each Riemann sheet This latter 
property also follows from eq. (|33|) since the integrand there is strictly periodic; this is a direct 



consequence of the sum over u n . As noted in |28[, a perturbative expansion in powers of A 



would display a spurious breaking of gauge invariance at finite temperature. This is because 
the perturbative condition, A small, is not preserved by large gauge transformations. This 
simply means that a truncated series expansion of a periodic function is not itself a periodic 
function in general. 

The behavior under parity is better exposed by introducing the functions <p n of Ap- 
pendix [B[ In the present case 

Q' a (uJ,0) = <j) (u) - au, 

0o M = log(2cosh(u;)) . (37) 

The parity transformation corresponds here to uj — > —uj. In agreement with our previous 
considerations, the term that breaks parity, — au, is a polynomial. Also, a parity transfor- 
mation is compensated by a change a — > —a, in agreement with eq. (0). The effective 
action can then be rewritten as 



W a (D) = -tr log(2 cosh(^(m + rjA ))) + rja^tv(A ) + aN f ^m 



(38) 



The first term coincides with the result in p8| after subtraction of a Ao-independent term 
(namely, the same expression with Aq — 0). It preserves parity but may break gauge 
invariance since it changes by integer multiples of in. This is compensated by the second 
term which is odd under pseudo-parity and breaks both gauge and parity symmetries. It is a 
local polynomial that introduces a parity anomaly. It cannot be removed without breaking 
gauge invariance. The last term is also a local polynomial that breaks parity but it can be 
removed preserving gauge invariance. This term illustrates that the ^-function prescription 
may yield unessential pseudo-parity even contributions to the anomaly; the same holds for 
chiral symmetry in even dimensions. The (essential) parity anomaly is thus 

A P = T]a J dx tr(A ) = rjaWcs , (39) 

where Wos is the 0+1-dimensional version of the Chern-Simons action. This checks that 
the essential parity anomaly is a local polynomial and temperature independent. When the 
gauge field configurations are traceless, Wcs vanishes identically. This is consistent with 
the fact that in this case the gauge group must be a subgroup of SU(iVj) which is simply 
connected and does not have large gauge transformations in 0+1 dimensions. (The gauge 
group itself may have large gauge transformations but the effective action must be strictly 
gauge invariant since one could always deform the configuration within SU(A^).) 
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The effective action can be written as the sum of components with well defined pseudo- 
parity, namely, W^(—A ) = ±W^(A ). After removing the term <jNf{3m/2, the pseudo- 
parity even component is gauge and parity invariant. The odd component can be identified 
as the anomalous part of the effective action at finite temperature in 0+1 dimensions since 
it contains all the anomalies 



^ anom (D) = W~(P) = || dxatrfa) , (40) 



where 



$ CT = |o;(^(m-A ),0)-p. P .c. 

= I [log (e-^-Ao) + i) - bg (e-^ m+A ^ + l)] . (41) 
P 

In this expression, p.p.c. stands for pseudo-parity conjugate, A — > —A . Besides, the branch 
cut of the logarithm is to be taken along the negative real axis to obtain the manifolds Ai a 
for each a. Alternatively W / an0 m(D) can also be written as 



^ anom (D) = 77(|Wcs + r odd ) (42) 



with 



odd 



-tr 



tanh-ftanh^W^' 



(43) 



This same function appears also in exact 2+1-dimensional results [[30[ (cf. subsection |IV Q ). 

The anomalous action W / anom (D) resulting from the calculation is to be compared with 
the naive temperature independent action 

l v (a-e(m))W cs , (44) 

which also saturates the parity and gauge variations of the full effective action but displays 
a singularity at m = which is spurious at finite temperature. In fact this is the zero 
temperature limit of W anom . On the other hand, the perturbative result follows from using 

4(7 

^-^r^J A + O(At). (45) 

Substituting in W anom yields 

WWD) = ^vU- tanh(^P)) Wcs + O(A 3 ) . (46) 

This is the usual perturbative result which violates gauge invariance under homotopically 
non trivial gauge transformations. 

It is noteworthy that in 0+1 dimensions the exact result can also be written in closed 
form without gauge fixing 

W a (D) = -trlog(e- a/3m fi r " T + 1) , ft = Te~tf Mxo)dai0 . (47) 
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IV. SUMMARY OF THE MAIN RESULTS AND DISCUSSION 



In this section we will summarize the results for the 2+1-dimensional case. The details 
of the calculation are given in next sections. The calculation proceeds by: (i) Expanding 
the expression in eq. (57) to select terms with at most two spatial covariant derivatives, (ii) 
Taking the Dirac trace. To make this step simple, we assume the scalar field M(x) to be 
a constant c-number mass m. (iii) Moving the operators Di to the right. This generates 
gauge covariant objects of the form [Di, X]. (iv) Carrying out the k integration. This kills 
all remaining gauge non invariant objects tr(0\XDi\0). This step requires to diagonalize Dq 
in the Hilbert space of flavor and time (with x fixed) and can be done in an Ao-stationary 
gauge, (v) Summing over fermionic frequencies u n . And (vi) integrating over z in closed 
form whenever possible. 

Wq(D) denotes the zeroth order term in the expansion in powers of D u W£(D) denote 
the pseudo-parity even and odd second order terms, respectively. The dependence on a and 
(3 will not be displayed explicitly. 



Wo(D) 



1 

47T 



(Txtl 



Ij m0i(|(m-A))) - (|j M^( m 



A))) 



-am 



+ p.p.c. 



(48) 

The functions (f> n (uj) are of the form <9~( n+1 ) tanh(u;), that is, they are primitives of tanh(a;), 
with appropriate integration constants so that they are even or odd functions of u. Their 
proper definition is given in Appendix [B| The function 4>q{(jj) has already appeared in the 
0+1-dimensional case. The notation p.p.c. stands for pseudo-parity conjugate, A — > —A , 
and after taking / <i 3 a;tr it coincides with the complex conjugate. 

2 



+ 




4)))A^A-|</>o(^(m 



dt tanh(|(m - t - A ))Ai l-V + — 



1 



2(m 



1 1 m\ 
-f p.p.c. (49) 



t)+V 

The symbol V denotes the covariant derivative T> (X) = [D ,X] and p.p.c. now includes 
V — > —Vq. In addition, we have introduced the new field Ai(x) which is defined as any 
solution (in the space of matrix valued functions, i.e., multiplicative operators in x and Xq 
spaces) of the equation 

V Ei = VIA. (50) 

Where Ei(x) = F 0i = [Do, Di] is the electric field. Essentially, Ai = Vq Ei. However, since 
the operator Vq is singular, V^Ei either does not exist or is not unique. Because Ai is not 
unique, we will have to check that the final result is independent of the particular solution 
taken. This check is done below. Clearly, the space of solutions of this equation is gauge 
and parity invariant. 



W 2 -(D) = £| d^eijtx 



AiAj)$ a - aAiV Aj 



+2m 



tanh(f(m-t + yl )) tanh(f(m - 1 - A )) 



2(m -t)+V 



2(m -t) -V 



(51) 
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The quantity $ CT was introduced in eq. (§!]) and Fij = [D i} D 



As already noted at the end of subsection |II Q , a further expansion can be done in powers 
of V . This expansion can be regarded as one in powers of = [D^, ], since Di appears 
only in commutators. Such an expansion preserves two essential properties of the original 
expansion in powers of Df. first, it does not break gauge or parity symmetries and second, 
higher order terms are increasingly ultraviolet convergent. 

Clearly, Wq is already independent of V . The expansion of W 2 is less obvious and 
is done in their respective sections. For the pseudo-parity even component we find (cf. 
subsection |VTT]) 



W^tleadinglD) 



1 



247T 



d 3 xtr 



tanh(|(m- A )) (3 



1 



2m 



2 4 cosh 



m 



A))) 



£? + p.p.c. (52) 



On the other hand, for the pseudo-parity odd sector, the leading order in V is (cf. sec- 



tion VII): 



W anom (D) 



8^ 



2 Fij A.iA.j ) $ c 



-~(tanh(^(m- A )) +tanh(^(m + A )) -2ajA i V A j 



(53) 



This component of the action differs from W% by terms which are ultraviolet convergent, 
gauge and parity invariant and free of multivaluation, thus it is natural to identify it with 
the purely anomalous component of the effective action: it is the ultraviolet logarithmically 
divergent part of the effective action that contains all the essential anomaly. 

In the remainder of this section we will discuss the properties of Wanom(D) since this 
term has attracted more attention in the literature. The discussion of Wq and will be 
made in their respective sections. Nevertheless, it can be pointed out that the pseudo-parity 
even components are strictly gauge invariant, not just modulo 2iri. This was to be expected 
since the only allowed ambiguity in the ^-function regularized action is 2nin which is purely 
imaginary and W + is real. Also they are parity invariant, except for the term — |crm 3 in Wq 
which is is removable by counterterms. On the other hand, W2 has the same parity and 
gauge transformations as Wa nom (D). No scale anomaly is present in the effective action. 

First of all, let us check that W / a nom(E ) ) does not depend on the particular solution 
taken for Ai(x). In order to show this, note first that if Aq(x) is not only stationary but 
also diagonal [A (x) , diA (x)] = 0, and a solution is simply Ai = A^. In an arbitrary 
A -stationary gauge, an explicit solution is Ai = Ai — UdiU^ 1 , where U(x) is any of the 
stationary gauge transformations bringing A (x) to diagonal form. Adding an arbitrary 
solution of T>o(X) = 0, the most general form of Ai is 



Ai(x) = Ai(x) - U(x)diU-\x) + Bi(x) , [A (x), Bi(x)] = 



(54) 



19 



Bi(x) is time-independent and commutes with Aq(x) at each point and is otherwise arbi- 
trary.^ An explicit calculation shows that, as expected, the field Ai transforms covariantly 
under stationary and discrete gauge transformations, modulo redefinitions of Bi. It is in- 
variant under parity transformations. Using antisymmetry of and the cyclic property of 
the trace, it easy to see that W / an om(D) is uniquely defined, that is, the dependence on Bi(x) 
cancels in eq. (|53|). We remark that Ai depends both on x and xq, in general. 



A. Parity variation 



Let us study the properties of W^anon^D) under parity transformations. As expected, 
replacing D by D p and o by —a leaves the formula unchanged thus eq. (|22|) is verified. 
Therefore, the anomaly can be obtained from the relation Ap(D) = W a (D) — W_ CT (D). 
Using the property 



- &- a = 4<t A , on M 
and substituting in Wanom(D), the parity anomaly is 



Ap(D) = t/o"^— J d 3 xeijtr 



(F i:j - 2A l A j )A - AiV A 3 



(55) 



(56) 



This expression can be simplified by using the following replacement (valid inside / d 3 x tr) 

eijAiVoAj = eijiAdoAj - 2A i A j A ) . (57) 

This follows from noting that AidoAj can be replaced by AidoAj since Ai and Ai differ only 
by a stationary field (cf. eq. (|54])) and integration by parts. Then, the terms containing Ai 
cancel and the final expression for the parity anomaly becomes 



A P (D) = J (PxeijtriAoFij - A^A^) = 7]aW C s(D)- 



(58) 



Wqs being the Chern-Simons action. This checks that the parity anomaly has a temperature 
and mass independent form and is a local polynomial, as expected from general arguments 
given in subsection [11 . 



B. Gauge variation 

Any correct calculation must be consistent with gauge invariance, however, since we 
have partially fixed the gauge, we can only check invariance under the restricted set of 
transformations considered in subsection |II C| , i.e., time- independent and discrete gauge 



6 This is the generic case. If two different eigenvalues of Aq(x) differ by 2irin/[3 in a sufficiently 
large region of the space-time, T>q(X) = will admit time-dependent and not Ao-commuting 
solutions as well. 
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transformations. In other words, it must be verified that bringing the gauge configuration 
to two different A -stationary gauges gives the same result modulo 2iri. As has been noted 
in subsection [11 Q for some gauge configurations all allowed gauge transformations (i.e., 
consistent with the condition doA Q = 0) are topologically small and on the other hand, going 
to a Ao-stationary gauge may require large transformations. Thus our calculation does not 
constitute a proof that the effective action is invariant under all gauge transformations. This 
fact has already been established as a consequence of using the ^-function renormalization 
prescription and the fact that gauge transformations do not change the spectrum of the 
Dirac operator. 

In order to study the gauge variation of W / anom (D), note that all the quantities there 
transform covariantly under time-independent gauge transformations. In particular, A{ is 
defined by a covariant equation. Regarding discrete gauge transformations, A is the only 
quantity which transforms inhomogenously, namely A^(x) = A (x) + A(ac). The terms 
with tanh in Wa, nom (D) depend on exp(/3A ) and this quantity is unchanged by discrete 
gauge transformations, cf. eq (|13|), hence the last term in W / anom (D) is invariant^. The 
only variation may come through $ CT . This quantity also depends on exp(/3A ), thus if 
the argument of the logarithm is chosen always on the same Riemann sheet, i.e., on the 
cut manifold M a , this quantity is trivially invariant. On the other hand, on the Riemann 
surface Ai one can consider a path connecting D and D^, for instance the path A = A + tA 
with < t < 1 and keeping m fixed. Let us consider first the generic non-degenerated case 
discussed in subsection |II C] for which eqs. (|T5|) apply, then a simple analysis shows that 

$^ - $ CT = 4a0(-(xm)A, on M . (59) 

This implies that the corresponding gauge variation of W / anom (D) is of the form (a — e{m)) 
times something independent of a. (Here 6 and e stand for the step and sign functions, 
respectively.) Because it has already been established in the previous subsection that the 
cr-dependent part of the action is ^rjaWcs, it follows that W / anom (D) must have precisely the 
same gauge variation as \rj(a — e(m))Wcs, that is 

W„{D U ) - iy anom (D) = W &(-am)(Wcs(V U ) ~ Wcs(D)) , on M . (60) 

Observe that when am > the formula predicts strict gauge invariance. This is correct 
since in this case the branch cut T and the spectrum do not intersect and so there is no 
flux of eigenvalues through the branch cut when going from D to D u . As noted previously, 
the variation of the Chern-Simons action is 2ni times an integer which depends on the 
topological numbers of the gauge transformation and the gauge field configuration. Then, 
as expected, both Wanom an d Ap are gauge invariant modulo 2ni. 

For degenerated configurations (defined in subsection [II C|) , A no longer needs to commute 



with Aq, however, by inspection of § a and unitarity of exp((3Ao), it follows that there will 
be no jump in the logarithm unless am < and in this case the jump is odd in a, therefore 
the same factor aO(-am) is obtained and the same argument as before applies. 



7 It is interesting to note that the two tanh terms in Wanom can be combined into 



sinh(/3m) 



which is the same combination appearing in [18| replacing Aq by the chemical potential 



cosh(/3m)+cosh(/3j4o) 
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C. Abelian-like reductions 



As a check of our formula, let us assume that A is also time-independent and commutes 
everywhere with Aq. In this case, the field Ai has (or can be chosen to have) the same 
properties. As a consequence all terms with Ai vanish due to antisymmetry and our result 
for WanomCD) reduces to 



W / an0 m(D) = ^ / ^Xe^tl 



(61) 



This expression holds also for W 2 ■ Separating the even and odd components under m 
—m, this can be rewritten as 



,(D) = vi^Wcs + r odd ) 



where 



r o dd — — — 



i 

47T 



(62) 



(63) 



is the expression introduced in Ref. [30|. As proven there (see also pp| , |31"| , |3"3"| ) , when A 
is in addition ^-independent, there are no higher order corrections and the right-hand side 



of eq. (plD gives the full pseudo-parity odd component of the effective action. (Recall that 
according to our definitions, r o dd is even under parity since m and A Q are both odd. r o dd is 
an odd function of m.) 



D. Zero temperature limit 

Using the zero temperature limit of $ CT , 

$ CT = 4a0(-am)A o , (T = 0) , 



(64) 



and the replacement in eq. (0), it is straightforward to derive the zero temperature limit 
of our formula, namely 



W anom (D) = \n{°- eH) ^cs , (T = 0) 



(65) 



This is the standard result at zero temperature (see e.g. 22 J33|). 



E. Massless fermions 

We may also consider massless fermions. In this case the terms with tanh cancel in 
Wanom(D)- Also, $ CT becomes an odd function of a on M.. Since all terms in the action are 



odd under a — > —a, it follows from subsection [IV A| that 



H / "anom(D) = -TjaWcs , for m = and on M . (66) 
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This formula holds for W2 too. This refers to the formula after analytical extension, i.e., on 
M.. Keeping track of the Riemann sheet in the logarithm in <fe a adds a term inn, so that the 
trivial gauge invariance of the ^-function prescription is preserved on M. a . The full result 
is an application of the Atiyah-Singer index theorem fl~4] . |29|j . Because in the massless case 



parity and pseudo-parity transformations coincide, W (D) equals \Ap{D). Thus, eq. (p6| 



is in fact exact to all orders |14j| . 



F. Perturbative result 



Next, let us show that our result for the pseudo-parity odd part reproduces the results 
obtained using perturbation theory at lowest orders. As noted in the Introduction, the latter 
approach yields a renormalization factor in front of the Chern-Simons action which is not 
quantized, and so breaks gauge invariance under large gauge transformations. To obtain the 
perturbative result we should retain terms of zeroth or first order in Aq in W an0 m(D)- Use 
of the replacement in eq. ([ 



i(D) = ~~ J (Pxeijtr 



and the expansion of $ CT in eq. fl45|) , gives 

4(7 



e <r/3m _)_ l v 2 



If- 



AiAj)A 



f3m s 



+ ^tanh( ^ 



a (AidoA, - 2AiAjA c 



0{Al) . 



(67) 



It trivial to check that the terms with Ai cancel, and the final result can be expressed as 



^ anom (D) = \vU~ tanh(^) J Wcs + 0(A 



0) > 



(68) 



which is the standard perturbative result [gj] . It has the same form as in the 0+1- dimensional 
case, eq. fliBT). 



G. Relation to the chiral case 



There are strong similarities with the situation of chiral anomalies in even dimensions. 
There the (consistent) chiral anomaly is defined as the variation of the effective action 
under chiral transformations. The chiral anomaly contains an essential part which can only 
be derived from a non polynomial action. As it is well-known, at zero temperature this is 
the gauged WZW action which again is pseudo-parity odd |§. As already noted, the chiral 
anomaly has also been shown to be temperature independent (see e.g. [0). This observation 
has lead to propose that the same gauged WZW action is the full anomalous action at finite 

The findings in 0+1 and 2+1 dimensions suggest that this is not the 
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temperature too 

case. Indeed, from the point of view of parity anomaly saturation and gauge invariance, 
the naive result given by eq. fl6"5D is entirely sufficient, however the correct result, eq. (|53|). 
has an explicit temperature dependence. It is worth noting that the naive formula has a 
singularity along the line m = which is spurious at finite temperature, since there are no 
infrared singularities, and is not present in the full result. If a similar situation takes place in 
the chiral case, the amplitude corresponding to anomalous processes (but not the anomaly 
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itself) will show a smooth temperature dependence. (By anomalous processes it is meant 
those processes driven by a logarithmically divergent pseudo-parity odd amplitude.) Such 
dependence has actually been found in Ref. |39|]. A study of the chiral case along the lines 



followed here has been carried out in |34 1 , with the result that the anomalous component of 



the effective action has indeed a non trivial temperature dependence. 



V. THE EFFECTIVE ACTION AT ZEROTH ORDER 

In this section we will compute the 2+1-dimensional effective action at lowest order in 
the derivative expansion. For simplicity we will use units (3 = 2. After setting Di to zero, 
eq. fl27|) yieldsQ 

Q s0 (D) = - / / — z s tr(0|— — ^ 10>. (69) 

S ' 0K ! J (2tt) 2 2 ^Jr 2ni v 1 loQ + IP + H 

Here Q = u n + D and fi = M(x) — z. As mentioned, we will take the scalar field M(x) = m. 
Then, Q, p and /i commute with each other and 70Q + TP can be brought to the numerator 

1 (Jt - 70Q - IP „ A _ . 2 n2,i,2 



(i = m-z, A = /i 2 -Q 2 + k 2 . (70) 



7oQ + IP + H A 
The Dirac trace can be computed immediately yielding 

^-/(g^SX^MOlflO). (71) 

Here the trace refers to flavor only. In order to perform the momentum integrals it is 
convenient to use integration by parts: 

—z s f(z) = — I —z s+1 f'(z) , (72) 
r 2tt2 V ; s + 1 Jt 2ni J K J ' v ; 

which holds provided that z s f(z) vanishes at infinite on T for sufficiently negative Re(s). 
Then 

The momentum integral can be done using 

f d2k k2n = T(n + l)T(N-n-l) N+1 
J (2tc) 2 (k 2 + M 2 ) N 4vrr(A) 1 ' { ' 



8 The subindex in f2 Si o(D) stands for zeroth order. We will no longer make explicit the dependence 

on r. 
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and it yields 



fi s „(D) = --E / — ^ +1 tr(0| (-? ^ IQ) . (75) 



The sum over frequencies can be done straightforwardly by using eq. fl32|) 

a^fp) = -±/ (_£_ _ _^_) tanh(/1 _ ^ + p . p . c . (76) 

Again, p.p.c. stands for pseudo-parity conjugate. In addition, we have used that A is time 
independent and so D can be replaced by A . 

Finally, the integration over z can be done using the functions f2r(u;, s) introduced in 
eq. (p4f) and Appendix O. Thus 



fi s ,o(D) = — / rf 3 xtr -fl r (m - A Q , s + 2) -fi r (m - A , s + 1) + p.p.c. (77) 

47TJ Vs + 2 s + 1 / 

To obtain the effective action it remains to apply — j-| s =o- As shown in Appendix M 



Q a (uj, n) = , n = 0, 1, 2, . . . , 



^,l) = 0iH-^(^ 2 -^(f) 2 ) , (78) 

Here the prime refers to derivative with respect to s and a = ±1 refers to the two possible 
inequivalent choices of T. The functions <f) n (uj) are defined for every integer n and satisfy 

<f>_ x (w) = tanh(w) , (79) 

0;H = n _!H, (so) 

n (-w) = (-l) n 0nH ■ (81) 

(See Appendix || for the fixing of the integration constants for non negative n.) The function 
<fio{uj) has already appeared in the 0+1- dimensional case and similar comments can be made 
here: for non negative integer n, the terms <f> n (u)) in Q' r (uJ,n) are those preserving parity 
and the breaking comes from the polynomial term which is odd under a — > —a. 

After taking the — ^\ s =o in eq. (|77|) and using eqs. (|78|), the effective action at zeroth 
order takes the form 

W (D) = y I d3xtT ( m <M m ~ A o) - 4>2{m - A ) - ^°™ 3 ) + P-P-C- (82) 

Recalling that m and Aq change sign under parity, it follows that this action preserves parity 
except for the last term am 3 , which is similar to the term am in 0+1 dimensions. Since this 
term is polynomial and strictly gauge invariant it can be removed from the action. 

Next we should check the gauge invariance of Wo(D). Under time-independent gauge 
transformations A changes by a similarity transformation and this leaves Wo(D) invariant 
due to the cyclic property of the trace. Under discrete gauge transformations the spectrum 
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of A changes by an integer multiple of in. The strict gauge invariance of f2 Sj o(D), and hence 
of Wo(D), on the manifold M. a (i.e., the complex plane cut along T) follows immediately 
from eqs. fl3"3|) and (|77| ) since tanh is periodic and so the functions Q-p(uj, s) are also periodic. 
The strict gauge invariance of W^o(D) also on the Riemann surface M. can be checked directly 
using the formula (cf. Appendix [TJ) 

4>n{u + y) - <Pn{u ~ y) = ine(m)— , n = 0,1,2,... (83) 

where we have used the notation m = Re(u>) (since in practice u = m =p Aq and Aq is anti- 
Hermitian), e(x) is the sign of x, and the difference refers to the straight path from u — y 
to uj + y on the complex plane. This identity can be used to show that the combination 
appearing in Wo(D), m(j)i{uj) — 02^) + c.c. (where c.c. stands for complex conjugate) is a 
periodic function. 



VI. THE PSEUDO-PARITY EVEN EFFECTIVE ACTION AT SECOND ORDER 



A. Dirac degrees of freedom 

Starting from eq. (pTI) , the second order takes the form (we use units (3 = 2) 



2 1 

I —p\ 10>. (84) 



Where Q and \x were introduced in the previous section. Using again eq. (|70|), the Dirac 
trace can be computed immediately. Since we are interested in the pseudo-parity even 
sector, W 2 + (D), we keep only terms without the Levi-Civita pseudo-tensor e^. The resulting 
expression is simplified using that kikj is equivalent to \k 2 5ij within the two-dimensional 
momentum integral. This yields 

^< D > = " / S? ? X tr<01 ^ * 2) X°4 

4^ fl «4A« W),,). (85) 

AAA A A A A A A/ 17 v ; 

Here the trace refers to flavor only. 



B. Space-time and flavor degrees of freedom 

Next, the spatial covariant derivatives are brought to the right producing the quantity 
Ei = [Q,Di] = F 0i . We have integrated by parts terms of the form [Dj, Ei] in order to 
produce a simpler expression. In this form the following expression is derived: 
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d 2 k 



E 



dz 
r 2m 



z>tr(0| 



+ A A 2 + A A 2 + A 2 A 2 + A 2 A 2 + A 3 A + 2 ^^^ 



V - A; 2 



i g q i q i 

A2 + A2 ' A2 + 2 A3 'A 




_ A A 2 + A 2 A 3 + A 2 A 3 + A 2 A 3 



g 



g 2 



A 2 A 3 



A 3 A 2 ' 



A 2 A 3 " 



|0). 



(86) 



Several remarks are in order here. First, the terms with explicit D t break gauge invariance 
and will be shown to cancel after momentum integration. That these terms break gauge 
invariance can be seen from Di\0) = Ai\0). The breaking is due to the state |0). The terms 
of the form EiEi are explicitly invariant under time-independent gauge transformations, 
since both Ej, and Q = D + uj n transform covariantly under such transformations. For 
these terms, the operators appearing in (0||0) are purely multiplicative in cc-space since all 
di operators appear inside commutators. Finally, to obtain a more compact expression, the 
cyclic property has been used in the gauge invariant terms. This is an important point of 
the present formalism and deserves further clarification. The gauge invariant terms contain 
the construction 



(X)=tr(0|X|0), 



(87) 



where X is purely multiplicative in a?-space. In this case, we can consider (X) for each x 
separately and use |0) to refer to the zero- frequency state in the Hilbert space of functions 
of Xq. The point is that ( ) is not a trace in rr -space and thus the cyclic property holds only 
in a restricted form. Nevertheless two cyclic properties can be used: 

Rule 1. IfX(Do) is a function of D only andY is a multiplicative operator in x-space, 



(XY) = (YX) . 



(88) 



Rule 2. If 'Xi(D ) and X 2 (D ) are functions of D only, andYi andY 2 are multiplicative 
operators both in x and xq the following substitution applies 

J2(X 1 (Q)Y 1 X 2 (Q)Y 2 ) = '£{X 2 (Q)Y 2 X 1 (Q)Y 1 ) , (89) 

n n 

provided that the sum overn is sufficiently convergent. (The sum refers to the n dependence 
in Q = D + uj n .) 

To see this, we exploit that A is stationary to introduce the basis of the flavor-time 
space (for given x) formed with eigenstates of Q, \a,£), namely 

A \a) = a a \a), (x \£) = e KeX ° (k, = mi) , 
\a,£) = \a) <g> \£) , 
Q\a,£) = (oj n + a a + K e )\a,£). (90) 
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Then, the first rule follows from 



(X(D )Y) = "£X{a a )(a,0\Y\a,0) = (YX(D )) . (91) 

a 

And the second rule comes from 

J2(Xx(Q)Y 1 X 2 (Q)Y 2 ) = E E E x iK + 0|H|/3, £)X 2 (u n + a p + k £ )(/3, £\Y 2 \a, 0) 

n n a,f3 e 

= EEE^K + fla - Kt)(a, -£\Y 2 \(3, 0)X 1 (u n + a p )(/3, 0|Ki|a, 
= T,T,T, X ^n + a a )(a,0\Y 2 \^£)X 1 (u n + a^ + K i )(/3,£\Y 1 \a,0) 
= y £(X 2 (Q)Y^X 1 {Q)Y 1 ) . (92) 

n 

The first equality uses that Yi }2 are multiplicative in a?-space. The second equality comes 
from exchanging a with (3 and £ by —I. The third equality uses that Yi >2 are multiplicative 
operators in xo-space and also n has been shifted to n + £. This last step means that the 
formula only needs to hold if Y^ n (Xi(Q)YiX 2 (Q)Y 2 ) remains unchanged under the replace- 
ment Q — > Q + u r (the same shift for all Q in the expression) which will be true if the sum 
over n is convergent. 

We remark that the present formalism actually only requires A (x) to be everywhere di- 
agonalizable but not necessarily ant i- Hermit ian or even normal. This is relevant if one wants 
to study the finite density case which requires introducing a chemical potential. Effectively 
this amounts to add a constant Hermitian term to A (see e.g. 1 37 . T 



C. Momentum degrees of freedom 

The momentum integration can be carried out immediately for the term with Df using 
eq. (|74|) and this term vanish identically. For the other non-gauge invariant terms of the 
form EiDi, we can replace by Ai since <9j|0) = 0. Thus all remaining operators in (0||0) 
are multiplicative in ai-space. 

In order to integrate over A;, we can again make use of the flavor-time basis \a,£) so 
that only c-numbers are involved. However, the integrand is a sum of terms of the form 
(XEiX'Ei) or (XEiX'Ai), where X, X' are functions of Q and this allows to use an equiv- 
alent and preferable method, namely, to use the label 1 to denote operators in the position 
X and the label 2 for operators in the position X '. More generally, operators X" in the 
position (XEiX'EiX") would also carry a label 1 due to the first cyclic property, eq. (j8SP . 
That is 

(XEiX'EjX") = (XtX'zXfEiEj) . (93) 
The second cyclic property, eq. (|8~9"D , then implies 

(Z(Q U Q 2 )EiEj) = (Z(Q 2 , Qi)EjEi) (94) 
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provided that Z(x\,x 2 ) = Z{x\ + u> r ,X2 + u) r ). With this ordering convention Qi and Q2 
are commuting objects and the momentum integrals in eq. fl86|) can be performed straight- 
forwardly. All required integrals can be derived from eq. fl?3]) and 



d 2 k 1 1 log(M 2 /M 



(2tt) 2 (F + M 1 2 )(A; 2 + M|) 4tt Mj 2 - M 



2 



(95) 



A non trivial check of the calculation is that all non-gauge invariant terms of the form 
(X EiX' Ai) cancel at this step. The result can be written in the following form 

where 

The symmetry of H(x\, X2) under exchange of its arguments is a direct consequence of the 
cyclic property, eq. (pp. 



D. The frequency degree of freedom 

In order to sum over frequencies, we first apply integration by parts, eq. (f72|) , to transform 
the logarithmic term in H{x\,X2) into a rational function. The resulting rational functions 
can then be reduced to a sum of simple poles and summed over n using the identity in 
eq. fl32|). The relevant formula is 



[iJ2H(x 1 + u n ,x 2 + uj n ) ~ (98) 

n 

( 11 z 11 z 2 1 \ tanhfu — xr) 

- H H 7 T7 7 — — + x i 2 + P-p.c 

\2 s + 1 xi — X2 {s + 1) (s + 2) xi — X2 J 2/j, — X\ + X2 

The symbol ~ is used since this substitution only holds within J r dzz s due to the integration 
by parts. X 12 means the same expression but exchanging the labels 1 and 2. Finally, p.p.c. 
stands for xi j2 — > —X\,i- When this formula is used in eq. (|9"E|), x\ = D i and x 2 = D 02 
(where D i and D 2 stand for D at positions 1 and 2 respectively) and so p.p.c. yields the 
pseudo-parity conjugate. 

In order to integrate over z it is convenient to extract the terms which diverge as a 
polynomial for large z: 



li H(x 1 + u n , x 2 + u n ) ~ - 



z 1 ls + 1 1 m 

+ 



s + 2x!-x 2 2 s + 2 (s + l)(s + 2)xi-x 2 



(ls + 1, , s 2 m 2 \ 1 

+ -{Xi - X 2 ) + 7 -T7 rrm + 



2s + 2 v (s + l)(s + 2) (s + l)(s + 2)xi - x 2 ; 2/x - xi + x 2 _ 

x tanh(yU — xi) + Xi j2 + p.p.c. (99) 
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E. The proper time degree of freedom 



Next we have to integrate over z and apply —d/ds\ s= Q. The terms in eq. (PUT) which are 
of the form (az + b) tanh(u; — z) are immediately worked out using the function fl a (uj, s). 
For the terms of the form tanhfy — z)/(u' — z) we will need some identities. Let f(z) be a 
meromorphic function such that zf(z) — > as z —>■ aoo, and g(s) analytic at s — 0, then 



u(») L^m) = ' ( 10 °) 

/ s=0 



r 2iri 



(] ( g{s) [ ^z s f(z)) =g(0) I dtf(t). (101) 



ds \ Jr 2ni 



s=0 



The first equality follows from taking s = directly in the integrand, since everything is 
convergent. The second equality follows from previous one. In our case, the two terms 
in eq. (|99|) which are of the form tanh(a; — z)/(u' — z) which are even on the explicit m 
are convergent for large z after adding their pseudo-parity conjugate, so eq. ( |101| ) applies 
directly to them. The remaining term is not directly convergent. It is of the form 

which can be rewritten as 

d ( , . f dz „ ( tanhfcij — z) — a o \ \ , 

h = ~Ts L [ l-2z + —;))^- (103) 



The subtracted term is convergent and thus it vanishes due to eq. ( |101|) . The remainder can 
be computed explicitly 

/ 2 = S(0)|. (104) 
The effective action takes the form 



W 2 + (B) = -Ur<0| 

S7T 



1 /11m 

- S) ; (m - Dol ,i) + ( i -- 5 _ 15 -)^(„,-A ll , ll 



am [-°™ ,.(l /r . n v . m 2 \ tanh(m - 1 - D 01 ) 

+ — h / at -^01 — + 



A y uzy D 01 -D 02 J 2(m-t)-D 01 + D 



02 



E 2 

|0)+X 2 + p.p.c. (105) 



(D i - A)2 



,2 



This expression can be transformed by making some observations. First, due to the 
cyclic property, the terms Xl j2 just give a factor 2 to the action. Second, recalling that 
V (X) = [Dq,X], the cyclic properties imply 

An - A)2 = V 01 = -T?o2 • (106) 

For instance 
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((An - D 02 ) 2 XY) = (DqXY - 2D XD Y + XD 2 Y) 
= (X[D ,[D ,Y}}) 
= {XVlY) 



(107) 



Therefore, D 02 can be replaced everywhere by D m + V Q2 . Third, because A commutes with 
d and d \0) = (0\d = 0, the symbol -D i can De replaced by A i everywhere. Fourth, we 
make use of the fields Ai(x) which were defined in eq. (|50|). The relation V Q Ei = T)\Ai 
allows to make the replacements 



E E ■ = — A A- E E 



1 



-AiV Aj . 



(108) 



(The minus signs come from by part integration.) Finally, we can make use of the functions 
(j) n {ui) instead of VL'^uj.n) (see Appendix |B|). After these manipulations, the pseudo-parity 
even effective action at second order can be written as 



W+(D) = j~j dx * tx 



1 



~0x(m - A )Ai—Ai - o (m - A )Ai - + -— A 



1 m 



4 



+ 



1 



dttanh(m - t - A )Ai -V + — — 

V4 V I 2(m-t) +V 



2 Do 



p.p.c. (109) 



(Here we are no longer using an ordering prescription with labels 1 and 2; the position of 
the operators is that given literally by the formula.) 



F. Transformation properties of the result and limit cases 

First of all, let us show that H^2 + (D) does not contribute to the gauge or parity anomalies. 
Because the equation defining Ai (eq. (0)) is gauge and parity covariant, we can choose Ai 
to transform covariantly too. The gauge invariance under time-independent transformations 
is immediate since A Q , Ai and T> transform covariantly. Under discrete transformations, 
At and T> remain invariant, whereas A changes by integer multiples of m thus the integral 
term in W / 2 + (D) is invariant due to periodicity of the hyperbolic tangent. For the terms 
containing <p n the invariance under discrete gauge transformations follows from eqs. (|83[) 
through a cancellation among the terms with <ft , <pi and their p.p.c. Since the cancellation 
is non trivial it provides a check of the calculation. 

The parity invariance of the terms containing <p n in VT 2 + (D) follows immediately from the 



parity properties of these functions, eq. d8T|) . On the other hand, in the term containing the 
integral over t, a parity transformation can be seen to be equivalent to change o — > — a, in 
agreement with eq. (p^). Thus the parity violating contribution would come from extending 
the integral on t to the range — oo < t < oo. Taking into account the p.p.c. term, such 
integral vanishes because it is convergent and the integrand is an odd function of m — t. 

Next, let us check that W 2 + (D) only depends on Ei(x) and not on the particular solution 
taken for the field Ai(x). Recalling the definition of Ai (cf. eq. ((50)) this requires that AiAi 



should appear with at least two powers of the operator V . The identity V 02 = Q 2 — Qi 
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(eq. ( |106|) ) then requires in eq. (|96|) that the function H(x\, X2) should be of order [x\ — X2) 



as x\ — X2 — > 0. In fact this is the case, 



H(xi,x 2 ) 3fi 2 + x 2 
lim 7 V2 = a~T~2 • ( 110 ) 



In eq. ( |109|) the dependence through Ei is not manifest because we have used the symmetry 



between the labels 1 and 2 to simplify the expression. To check this property directly from 
the final expression for W 2 + {J^>) it is necessary to restore D01, -D02, add (dividing by two), 
eliminate Dq 2 in favor of T>o and then expand in powers of T>q. (This procedure is carried 
out in detail for W 2 ~ in section |VH] .) After doing this, it is indeed found that the necessary 
cancellations take place and only terms of order A{D^Ai survive. These cancellations among 
different terms are again a non trivial check of the result. Retaining the lowest order in the 
expansion on V just described yields 

, /T ^ N 1 [ ,0 rtanhfm — An) 1 1 9 , 

^ lcadmg (D) = - — / d*xtr \ -Bl - 7 —^ 7 - - - 1^ + p.p.c. (Ill) 



2m 4 cosh (m — An 



It is noteworthy that the terms in the expansion combine in such a way that the integral over 
t present in the original expression can be done in closed form by integration by parts. It 
has been checked that the same final expression is obtained starting directly from eq. (|110 ). 



In fact, the interest of W 2 leading goes beyond checking explicitly the independence of 
W 2 + on the spurious degrees of freedom contained in the variables A4. W 2 + t leading has been 
obtained by making first an expansion in powers of Di to second order and then keeping 
the leading order in an expansion in powers of T>n. As noted at the end of subsection [II G| , 
this double expansion in Di and T>q is the natural finite temperature version of the standard 
space-time gradient expansion at zero temperature (which in the absence of other non-gauge 
fields, is also equivalent to a 1/m expansion). No symmetry has been broken by the further 
expansion in powers of V and W^tieading i s manifestly gauge and parity invariant. On the 
other hand, expanding in powers of Aq, as implied by standard perturbation theory, would 
break invariance under discrete gauge transformations. 

In the particular case of Abelian and stationary fields (i.e., time-independent A) T> 
vanishes so W^t leading does in fact coincide with W 2 + ■ 

It is also of interest to examine the zero temperature limit of W 2 + (X>)- As shown in 
Appendix [B], 

lim (pn-Auu) = — t u n e(m) , n = 0,1,2,.... (112) 

(Recall that we are using units (3 = 2 so the zero temperature limit corresponds to large 
uj = m — Aq). After some non trivial cancellations, one finds 



D) = d 3 xtrEi —- + -Urn 2 + V 2 A dt - - 

4tt J 8\m\ 2 K 0J J\ m \ At 2 {At 2 -Vl) 



E iy (T = 0) . 

(113) 



This result depends only on Ei (and not Ai) and is free of parity or gauge anomalies. Full 
Lorentz invariance, which is recovered at zero temperature, is no supported since it has been 
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broken by making an expansion on the spatial gradient only. In order to compare with direct 
zero temperature calculations, we will only retain terms of dimension 4, or equivalently, the 
leading term in a 1/m expansion. This gives 

^ = -^/*^ +0 (^)- (r - 0) - (114) 

(Of course, the same result is obtained taking the zero temperature limit of W^ieadiagi^) ^ n 
eq. (|llip ). Then Lorentz invariance implies 

W+(D) = --L £ -^ [d*xtrF^ + o(±) , (T = 0). (115) 
487T m J p \mv 

(Up to constant terms coming from Wo(D).) It is an interesting check of the calculation 
that precisely the same result is obtained from a direct zero temperature calculation (using 
for instance the 1/m expansion coefficients given in 

VII. THE PSEUDO-PARITY ODD EFFECTIVE ACTION AT SECOND ORDER 

In this section we will compute the leading contribution to the imaginary part of the 
effective action, which comes from the second order in the gradient expansion. This leading 
term is particularly interesting since it contains the essential parity anomaly as well as the 
multivaluation under large gauge transformations. 

Many of the ideas needed to carry out the calculation have already appeared in the 
previous section, so they do not need to be repeated here. We start from eq. (recall 
that we use units (3 = 2). Use of eq. (f70D allows to compute the Dirac trace. Keeping just 
the terms with yields 

+(/i 2 _ fc2) i D Q D . _L _ ^| A Ini _ D 1 D 1 \ | Q) (116) 

VP J A A A A A A A A J AJ 1 1 v ; 

where the trace refers to flavor only. (This expression is analogous to that in eq. fl55] ) for 
the pseudo-parity even part.) 

Next, the spatial covariant derivatives are brought to the right. It has been found to be 
best to move Di first and then Dj to generate a smaller number of terms. This involves the 
quantities Ei and Fij = [Di, Dj], and also A{ appears due to Dj|0) = Aj|0). In this form the 
following expression is derived: 
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~~K E Ti E + - ±E%E - S±E9-E - %E%E 

A A 2 A A 3 A A 3 A 3 A A A 3 

+(? + 2k^E^E - 0- + »») JjsglS - Sh E % E 

Q E % E _ + % E %E + k *± 2 E% E 

A 2 A 3 A 2 A 3 A 3 A 2 A 2 A 3 

+A; A^A^ - A E & A ~ A 2 A 2 " A^A^J |0) " 

This expression is analogous to that in eq. fl86l) for Qf 2 (D) an d many of the remarks made 
there apply here too. A standard differential form notation has been adopted for the spatial 
indices, i.e., dx\ and dx 2 are anticommuting, d 2 xe^ = dxidxj, D = Didxi, A = Aidxi, 
E = Eidxi and F = D 2 . The operator inside (0||0) is multiplicative in as-space, thus the 
part of (0||0) corresponding to that space has been replaced by J d 2 x and d 2 x has been 
included in the differential forms, so |0) now refers to the xo-space only. The terms with 
explicit A break gauge invariance and will cancel later. Finally, to obtain a more compact 
expression, the two cyclic properties have been used. 

To carry out the momentum integration, we again make use of the method of labelling 
the Q operators with labels 1 or 2 according to their position. Note, however, that now 

(Z(Q 1 , Q 2 )E 2 ) = -(Z(Q 2 , Q 1 )E 2 ) , (118) 

the extra minus sign coming from E being a 1-form. Since all quantities are now commuting 
(or ant i- commuting) momentum integration is straightforward. A non trivial check of the 
calculation is that all non-gauge invariant terms (X\EX 2 A) cancel at this step. The result 
can be written in the following form 

where H(x\, x 2 ) is the same function as in the pseudo-parity even casef], eq. (|97|). Note that 
in this case E 2 is odd under exchange of the labels 1 and 2 (since E is a 1-form), so again 
the total expression is even under this exchange. 

In order to sum over frequencies we need a new identity for the term with F, namely^ 



9 This suggests that the second order pseudo-parity even and odd components are closely related 
and perhaps connected through some kind of analytical extrapolation. Such relations exist in 1+1 
dimensions and allow to exactly compute the corresponding Weyl determinant |4(J . 

10 The series in the left-hand side is not absolutely convergent, nevertheless it will be so upon 
subtraction of a suitable ^-independent term, as in eq. ( j30|) . Thus this equality holds inside J r dzz s . 
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E 



X + UJ T: 



(x + u n y 



— tanh(/i — x) — p.p.c. . 



(120) 



In the term with H{Q\ 1 Q2) we again use integration by parts to eliminate the logarithm. 
This yields 



H(xi + u n , x 2 + u) n ) 



+ -x 



1\ tanh(/i — x\ 



S + 1 X\ — X2 2J2jJ, — Xi+X2 



+ Xi >2 + p.p.c. . 



;i21) 



Note that this result multiplied by \x = m — z needs to coincide with that in eq. (j98l) only 
inside / r dzz s , since integration by parts has been used in both cases. In analogy to eq. ([ 
it is convenient to rewrite this expression as 



Y^,H{xi + u n ,x 2 + 



UJr, 



1 1 



+ 



+ 



in 



2 s + 1 X\ — X2 \2 s + 1 s + 1 X\ — X2/ 2/i — X\ + X2 
x tanh(/i — x\) + X 1>2 + p.p.c. (122) 

After carrying out the ^-integration, we arrive at a formula analogous to eq. ( 105|) , namely 



W 2 "(D) = £tr(0| 



ln' a (m - D , 0)F + J f J^-0^(m - D 01 , 0) 



+ 



(7 



m 



4 P 02 



tit 



2\2V 02 
tanh(m — t — Dqi) 
2(m -t)+ V 02 



E 2 



n 



02 



|0) - p.p.c. 



(123) 



The contribution refers to everything inside the parenthesis and it gives a factor 2 due 
to the cyclic property. The final form is obtained by further replacing Dqi by Aoi, E 2 by 
-Vl 2 A 2 and (0|0) by / dx . This gives 



W 2 -(D) 



8^ 



+2m 



dxntr 



(F - A 2 )$ a - oAVqA 

7 tanh(m — t + A ) tanh(m — t — A ) 



dtA 



2(m -t)+T> 



2{m-t) -V 



A 



(124) 



where 3> CT = Q' a (m — A , 0) — p.p.c. was introduced in eq. (fUD. Note that d 2 x is included in 
the differential forms. 

In order to obtain the leading term in a further expansion in powers of T>q (the analogous 
of W^ieadingC-D)) we go back to the symmetrized formula, eq. ( |123|) , then express the operator 
Dq2 in X 12 as D 01 + V 02 and carry out an expansion in powers of V 02 keeping the leading 
order. At the end D m can be replaced by A m . This needs only be done in the term containing 
the integral over t: 
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m 

m 

m 
1 



dt 

i 

dt 
dt 



tanh(m — t — D i) 
2{m-t)+V 02 ' 

tanh(m — t — D i) 
2{m-t)+V 02 ' 
tanh(m — t — Dqi) 
2{m-t) 2 



-^1,2 

tanh(m 



t — D 



01 



V, 



in, 



+ 



2(m - i) 
1 



n 



02 



cosh (m — t — D 01 ) 2(m - t) 



V 2 + O(P 2 2 ) 



tanh(m - D 01 )V 02 + 0(V[ 



02 ' 



(125) 



So finally the result is 

W„(V) = l ir I dx Q tr (F - A 2 )$ c 

07T J 

+ -^tanh(r 



Aq) + tanh(m + A ) - 2a)AV A 



(126) 



Let us finish by discussing the properties of W 2 ~ paralleling the discussion of W anom in 
section |V|. Using the same arguments as for Wanom it follows that W 2 ~ is well defined, i.e., it 
does not depend on the choice of Ai(x). This implies that it is possible to express all terms 
not contained in W^ncm in terms of the electric field. 

The term of W 2 with the integral over t is actually independent of a since using +oo 
or — oo as the upper limit of the integral gives the same result. This follows because the 
integrand is an odd function of m — t. Therefore it is parity invariant and W 2 has the same 
parity anomaly as Wanom- Likewise, the difference between W 2 and PUanom is strictly gauge 
invariant. 

Eq. (|6"T| ) for Abelian-like configurations also holds for W 2 , since the terms without F 
cancel. In the massless case, the term with an integral over t cancels and so eq. ( |66j) still 
holds for W 2 . 

The zero temperature limit of W 2 is given by 



W 2 -(D) = r]aQ(-am)W C s - ^ / dx tx f°° dtE 

07T J J\m\ 



rnDc- 



t 2 (4t 2 



-E, (T = 0). (127) 



The first term was already present in W anom . Since the second term is of order 1/m 2 in an 
inverse mass expansion, it is ultraviolet finite. 

The zero temperature limit of W^~(D) is not directly Lorentz invariant. Restoring this 
invariance requires including higher order terms in the gradient expansion. Nevertheless, 
keeping the first non trivial 1/m contribution in eq. ( J127|) , the only Lorentz invariant com- 
pletion is 



W~(D) = W Q(-am)W cs - £f^r / (^V^W^) + O (J^j , (T = 0) 



128) 



Precisely the same result is obtained through a direct zero temperature calculation using 
the formalism of Ref. 1 1231 1 . This is another non trivial test of the calculation. 
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APPENDIX A: THE ^-STATIONARY GAUGES 

Let the space-time have the topology M^+i = M d x S 1 , d being the space dimension. All 
functions involved are assumed to be continuous on M d +i (and thus, in particular, periodic 
in the time coordinate) unless otherwise stated. Let A^(x) be a gauge configuration on 
M d+1 . We want to show that there exists another configuration B^(x) in the A -stationary 
gauge (i.e., d B Q = 0) related to by a gauge transformation. To do this, let us construct 
the following time-independent field 

Q(x) = Texp A (x,x )dx S j . (Al) 

(The symbol T stands for time ordered product.) As is well-known, Q transforms covariantly 
(in fact tr(Q) is just the Polyakov loop) 

Q u (x) = U-\x, 0)Q(x)U(x, 0) . (A2) 

Thus, Q is gauge invariant modulo stationary gauge transformations. In view of this, we 
demand that should yield the same field Vt as A^, that is, 

Q(x) = exp (-PBq(x)) . (A3) 

That the function B (x) above exists (with continuity on M d ) follows from Q being an 
homotopically trivial application from into G. In fact, 

U x {x,t) =Texp^- j* A Q (x,t')dt'^ , < t < (3 (A4) 

defines one such deformation of Q to the identity element of G with t as interpolating pa- 
rameter. It only remains to show that there actually exists a gauge transformation bringing 
A^ to Bfj,. This is achieved with 

U(x,x ) = Ui(x,x ) exp (x B (x)) . (A5) 

The first factor, U 1: brings A to zero but breaks periodicity in Euclidean time. The second 
factor e x ° B ° yields Aq = B and reestablishes the continuity of U on M d+1 . Therefore, a 
Ao-stationary gauge always exists. 

Next, let us determine the remaining gauge freedom within an Ao-stationary gauge. Let 
A^ and B^ be two A -stationary configurations and let U be a gauge transformation bringing 



A^ to B^. Then 



B {x) = U- 1 (x)d Q U(x) + U- 1 (x)A Q (x)U(x) . (A6) 
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In terms of the auxiliary variable V(x) = exp(xoA (x))U (x) , the previous equation becomes 

B (x) = V- 1 (x)d V(x). (A7) 

Its general solution is 

V{x) = U Q (x) exp(x Q B Q (x)) , (A8) 

where U is an arbitrary time-independent gauge transformation. In terms of U (x) 

U(x) = exp(-x A (x))U (x) exp(x B (x)) . (A9) 

This is the most general "gauge transformation" bringing into a A -stationary gauge, 
and depends on the two arbitrary time-independent fields B and U . However, U will not be 
a true gauge transformation unless it is also continuous on M d+1 . To impose this restriction, 
let us change variables from B (x) to A (a?) by means of 

B (x) = U 1 (x)(A (x) + A(x))U (x) . (A10) 
This allows to write U (x) more conveniently as 

U(x) = exp(-x A (x)) exp(x (A (x) + A(x)))U (x) . (All) 
The continuity of U(x) follows from requiring U(x,0) = U(x,f3), which is equivalent to 

exp(/M (aO) = exp(P(A (x) + A(x))) . (A12) 

Finally, let us show that for simply connected gauge groups, it is not always possible 
to bring any given gauge field configuration to the Ao-stationary gauge using only small 
transformations. To show this, let us assume that the initial configuration A^ is already 
Ao-stationary. Let us perform a large gauge transformation on it, yielding the configuration 
A" , and then bring this configuration again to an A -stationary gauge to finally obtain a 
configuration A' . Without loss of generality it can be assumed that A^ and A'^ are related 
by a discrete transformation of Aq(x) (i.e., with no time-independent part). We want to 
show that the gauge transformation bringing A"^ to A'^ cannot always be chosen to be small. 
Indeed, if it were small, A^ and A'^ would be related by a large gauge transformation which 
is also a discrete transformation for A (x). However, this not always possible if the gauge 
group is simply connected. This can be seen by taking A (x) to be a constant diagonal 
matrix such that the spectrum of exp(f3A (x)) is not degenerated. Then, the allowed A (a;) 
for this configuration will be also diagonal and constant and thus e x ° A will necessarily be 
small, that is, no large discrete transformation exists for such A (x). On the other hand, 
in the Abelian U(l) case, it is always possible to deform A (x) to a time-independent one 
preserving the Polyakov loop and no large transformations are required. 

APPENDIX B: AUXILIARY FUNCTIONS 

The functions D,r(u, s) are defined through the formula 
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f dz 

Q v (uj, s) = - —z s tanh(u; - z) . (Bl) 
Jr 2m 

The integral is defined for Re(s) < —1 and extends to a meromorphic function of s with 
a simple pole at s = — 1. The integration path T follows the ray 9 starting from infinity, 
encircles zero clockwise and goes back to infinity along the ray 9 — 2n on the z-complex 
plane. Because the function tanh(z) has simple poles at z = uj n = m{n + |), the integral 
is not defined when uj + uj n lies on V for some n. As a consequence, the complex plane 
corresponding to the variable uj is cut along rays with angle 9 stemming from the points 
uj = uj n which are branching points of the function s) (unless s is an integer). Let 

us denote the a;-complex plane so cut by Cg, and the corresponding Riemann surface by C. 
This latter manifold is common to all values of 9. 

The angle 9 can take any real value excepting 9 = ±7r/2 (mod. 2ir) (i.e., T on the 
imaginary axis). The function z s is defined taking arg(z) in the range ]9 — 2m9[. Corre- 
spondingly, we will use the notation (z)g, arg e (z) and flg(u, s). Applying Cauchy's theorem, 
the function can then be written as 

n e (uj,s) = J2(u-Lo n y e . (B2) 

n€Z 

If 9 is shifted to 9' but without crossing the limits ±7r/2 (mod. 2tt), only a finite number 
of terms in the series get modified and the two functions Qg, Qgr are related by analytical 
continuation, that is, they coincide on the C. On the other hand, from the definition it 
follows immediately that 

£Wn(^, s) = e 27Tins n e (uj, s), (ne Z) (B3) 

As will be shown below, fip(u;, k) = for k = 0, 1, 2, . . ., so Q'g^nni 00 ' ^) = ^0(^7 ^) (where 
the prime refers to the derivation on s). Thus, all angles differing by an integer multiple of 
27r give the same effective action and are equivalent. It follows that all practical cases are 
covered by taking 9 = tt and 9 = 2n, which correspond to a — 1 and a = — 1 respectively. 

The tu-complex planes C CT are cut along rays which are parallel to the real axis, start at 
uj = u n and go to the right for a — +1 and to the left for a = — 1. To describe the functions 
fl a (u>, s) on the Riemann surface C, it will prove convenient to introduce the a;-complex plane 
cut in another way, which will be denoted by C p and is defined as follows. Starting from any 
of C a planes, rotate upwards the rays which are above the real axis so that these rays lie now 
on the positive imaginary axis, and likewise rotate downwards the rays below the real axis 
to put them on the negative imaginary axis. Since the branching points which are nearest 
to the origin are uj = ±ujq = ±27r/2, the functions Q a (uj,s) extended to C p are analytical 
for all uj except when uj is purely imaginary and above ujq or below — ujq. The plane C p has 
the property of supporting parity transformations, which correspond to uj —>■ —uj (since in 
practice uj = m ± A Q ) and also discrete gauge transformations, uj — > uj + in, where the two 
points are connected by a straight line. Besides, for each value of cr, the function fl a (u), s) 
takes the same value on C CT and C p when uj is in the half-plane cxRe(co') > 0. 

The functions fi^a;, s) can be directly related to the Hurwitz function ((z,q) [41| (p. 
1073 and ff.): 

00 

C(z,q) = Y,(Q + n)-\ (B4) 

n=0 
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where the subindex ir means | arg(g + n)\ < tt. This function is analytical (as a function of 
q) except on the negative real axis. The series representation of Q a , eq. ( |B2|) , allows to write 



Sl+fa s) = (-tn)X (s, + MX (s, 1 + 

V 2 m) V 2 if/ 

fi_(w, s) = e^ s {- m )X (s, ]:--) + MX (-8, \ + -) , (B5) 

V I in J \ A ITT/ 



where all arguments are to be taken on ] — tt, 7r[ and the equalities hold on C p . An immediate 
consequence is that the functions Q + (u,s) and Q_(uj,s) coincide everywhere when s is an 
integer. Using the identities El| 



C (_ n;g ) = _^±lM ) Bn {\ - x) = (-l) n B n (x) , n = 0,l,2,... (B6) 



where B n (x) are the Bernoulli polynomials, it follows immediately from eq. (p5| ) that 

Q ± {u,n) = 0, n = 0,1,2,... (B7) 

and furthermore, 

n' + (0,n) = n'_{0,n), n = 0,2,4,... (B8) 

^(0,n) = -a(mr +1 ^±i^, n= 1,3,5,... (B9) 

n + 1 

where the prime refers to derivative with respect to s. Another useful property follows from 
the identity 

d u Q! a {u, s) = n a (oj, s - 1) + sQ' a (u 7 s - 1) , (BIO) 

which implies 

9X(w,fi+l) = (n+l)fl>,n), n = 0,1,2,... (Bll) 

The functions f2^.(ci;,n) for non negative integer n can be expressed in terms of simple 
integrals as follows 

Q' (oj, n) = — — f / z s tanhfu; — z) ) 

\ / s=n 

P — (TOO 

= / dtr(tanh(cj - t) - a) 
Jo 

= dt(u-t) n (timh(t)-(j), n = 0,1,2,... (B12) 

J (TOO 

The defining integral has been converted into a convergent one by using eq. ([JT]) and then 
eq. (|101| ) has been applied. This result refers to to in and in the integrals over t, the 
integration path is to be taken parallel to the real axis. (A similar method would also yield 
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the relations in eqs. (p7|JB8| , p9|) .) In particular, for n = the integral can be done in closed 
form 

fi>,0) = log(e- 2 ™ + l), (B13) 

which can be rewritten as 

n;(w,O) = 0o(w)-<Tu;, (B14) 

0o M = log(2 cosh(cu)) . (B15) 

Next, let us introduce the functions 4> n {oS) for arbitrary integer n. They are defined on 
C p by the relations 

= tanh(cu) , (B16) 

^H = 0n-iH, (B17) 

4> n {-u) = (-l) n n (u;) , (B18) 

0„(O) = i-fi' ± (O,n), n = 0,2,4,... (B19) 
n! 

The functions n for non negative n follow from eq. ( |B17| ), 4> n {uj) = n (O) + Jq dz <fi n -i(z), 

where the integration path lies on C p . For odd n, 0„(O) = due to parity, eq. (P18|) . For 
even n, n (O) is well defined from eq. ( |B19|) due to eq. ( p8|) . 

The relation between the functions <f) n (uj) and Q' ± (uj,n) can be established by starting 
from eq. ( |B14|) and recursively applying Jq duj. This yields 

— n' a (u,n) = <t> n {u) ~ crP n +i(v) , n = 0,1,2,..., on C p . (B20) 
n! 

Where P n (u) are polynomials of degree n which satisfy 

P (u) = l, (B21) 

P' n+1 {u>) = P n (u) , (B22) 

P n (-uj) = (-l) n P n (oo) , (B23) 
(in) n 1 

Pn(0) = 5n( 2 ) ' n = ' 1 ' 2 '--- ( B24 ) 



The lowest order relations are 



K 




o) = 




— au 


K 




i) = 




-( 






2) = 


2 (w) 


-J 



1 „, , „x /l 3 1 /ztt n 2 



2! ffV ' ' T y ' V 6 6 V 2 / 



(B25) 

) , (B26) 
w J . (B27) 



These relations make explicit the transformation properties of the functions f^(u;,n) 
under parity, u — > — cj. The polynomial part is responsible for an anomalous violation of 
parity. Also, it follows 
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tf (T (u,n) = (-l) n tf_ < ,(-u,n) (B28) 

that reflects eq. fl22"D . 

Next, let us discuss the periodicity properties of these functions, which are related to 
the discrete gauge transformations. Since tanh(o;) is a periodic function with period in, so 
are the functions f2o-(^, s) on C CT . This reflects the trivial gauge invariance of the (^-function 
regularization due to the invariance of the spectrum of the Dirac operator. Let us consider 
now the periodicity on C p . Since Q a (u, s) takes the same values on C a and C p on the half- 
plane cxRe(co') > 0, it follows that this function is periodic on that half-plane of C p . For 
the other half-plane, aKe(uj) < 0, take u in the strip — | < Im(u;) < | (which is common 
to Co- and C p ). Then, the value of tt a (u) + in, s) computed by taking u + in in C p will be 
different from the value arrived at on and this difference comes solely from the term 
n = in eq. (|B"2|). For instance, for a = +1, the contribution of this term on C a would 
be [uj + in — ujo)^, whereas on C p the contribution is instead (u + in — uo) s _ n - Noting that 
fl a (uj + in, s) = D,o-(lj, s) on C CT and also that Q a (u, s) takes the same value on C CT and C p 
because of — | < Im(u;) < |, it is easily established that, for to + in G C p , 

n a (u + in, s) - n a (u, s) = ( e - u2nis + uj ) s a &(-(rm) . (B29) 

Here ujq = in/2 and m = Ke(uj). After analytical extension beyond the strip — | < Im(cj) < 
|, this relation holds on the whole plane C p . From here it follows the useful relation 

fl' a (uJ + ujq, n) — Vt' a (uj — ujq, n) = — cr2ni u n Q(— am) , uj G C p , jiGZ. (B30) 

Where 0(x) is the step function. Using that, on C p , 4> n (uj) is the component of Q! a (uj,n)/n\ 
which is even in a, we also find 

4> n (uj + uj ) - <p n {u - uj ) = in —e{m) , uj G C p , n = 0,1,2,..., (B31) 

where e(x) stands for the sign of x. 

Finally, let us study the zero temperature limit properties of these functions. This 
corresponds to the large u limit. Again we start with uj in the strip — | < Im(c<j) < ~, then 
from eq. ( pi2|) it follows that 

~o;-,oc -2a -O(-am), lu G C p , n = 0,1,2,... (B32) 

n + 1 

By analytical extension, this relation holds on C p . Moreover, for positive am, the zero limit 
has only exponentially decreasing corrections. Then, from eq. (|B20p it follows that 



n -iM ~u,-*oo £(ra)P„(u;) , wgC p , n = 0,1,2,... (B33) 
and in particular 

^(u) = ^e(m) ( 1 + O ( -)) , ojeC p , n = 0,1,2,... (B34) 
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